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ABSTRACT 

Aims. We study the effects of local environmental conditions affecting the diffuse interstellar band (DIB) carriers within the Upper 
Scorpius subgroup of the Sco OB2 association. The aim is to reveal how the still unidentified DIB carriers respond to different physical 
conditions prevailing in interstellar clouds, in order to shed light on the origin of the DIB carriers. 

Methods. We obtained optical spectra with FEROS on the ESO 1.52m telescope at La Silla, Chile, and measured the equivalent widths 
of five DIBs (at 5780, 5797, 6196, 6379, and 6613 A) as well as those of absorption lines of di-atomic molecules (CH, CH + , CN) and 
atoms (Ki, Cai) towards 89 targets in the direction of Upper Scorpius. We construct a simple radiative transfer and chemical network 
model of the diffuse interstellar medium (ISM) sheet in front of Upp Sco to infer the effective radiation field. 

Results. By measuring the DIB and molecular spectrum of diffuse clouds towards 89 sightlines in the Upper Scorpius region, we 
have obtained a valuable statistical dataset that provides information on the physical conditions that influence the band strengths 
of the DIBs. Both the interstellar radiation field strength, / uv , and the molecular hydrogen fraction, / H , , have been derived for 55 
sightlines probing the Upp Sco ISM. We discuss the relations between DIB strengths, CH and CH + line strengths, E(b-v>, Ajv> and 
/h,. The ratio between the 5780 and 5797 A DIBs reveals a (spatial) dependence on the local environment in terms of cloud density 
and exposure to the interstellar radiation field, reflecting the molecular nature of these DIB carriers. 

Key words. Astrochemistry - ISM - ISM: lines and bands - ISM: molecules - ISM: dust, extinction - ISM: clouds - ISM: individual 
objects: Upper Scorpius 



1. Introduction 

The diffuse interstellar medium contains compounds of uniden- 
tified origin that absorb in the UV-visual to near-infrared spec- 
tral range. More than 300 di fferent diffuse interstellar band s 
(DIBs) are currently identified dHerbidl995l:lHobbs et al.l2008l) . 
Many possible carriers have been proposed, ranging from 
grain impurities and exotic molecules to H2. In the past two 
decades the field has converged towards larger carbonaceous 
molecules, like the fullerenes and polycyclic aromatic hydro- 
carbons (PAHs), w hich have electronic transitions in the opti- 
cal (see for example Salama et al.lll999llRuiterkamp et al.ll2.QQ5t 
Kokkin & Schmidt! 120061 Zho u et al 1 120061 and Salama et al. 



In order to understand the chemical and physical proper- 
ties of the DIB carrier(s) it is important to study their be- 
haviour in different interstellar environments, both in our own 
galaxy and beyond . Stud i es of DIBs in the Magel l anic Clouds 



(lEhrenfreund et all 20021 Cox et al 
20061). M31 (ICordiner et al.l 12008 



2006L I2007t IWeltv et al.l 
IT 



- . affbl) an d beyond (e.g. 
Heckman & Lehnertl 120001: lYork et al.l 120061: ISollerman etatl 



20051: ICox & Patajl2008h illustrate that DIB carrier abundances 



201 lh . New diffuse bands have been detected in one line-of 



(per amount of dust and gas) can be similar to galactic values. 
However, these studies have also revealed systematic differences 
in these extragalactic environments. 

A large amount of published information is available re- 
garding DIBs in many sightlines probing the Galaxy (e.g. 



sight which appear to match with naphthalene and anthracene 
cations dlglesias-Grofh et al.ll2008l 1201 Oh and the weak 5450 A 
DIB is found to match with an absorption ban d arising from 
a hyd rocarbon plasma created in the laboratory (iLinnartz et al.l 
l2010l) . Linear-C3H2 has been put forward as a carrier of the 5450 
and 4881 A DIBs bv lMaier et al.l d201 lh. These assignmen ts are, 
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properties with respect to each other and to other diffuse ISM 
gas & dust tracers. Most of these studies focused on DIBs prob- 
ing various galactic environments, and provided average results 
for the Milky Way. Studies dealing with a particular region usu- 
ally only include a very limited number of sightlines. An excep- 
tion is the study of the Orion region bv lJenniskens et all (1 1994b 
which entails 22 line s of sight. Another multi-object study, by 
Ivan Loon et a II (120091) . used the distant globular cluster a> Cen to 
probe fluctuations of Can, Nai and the A/15780 and 5797 DIBs 
in the diffuse - low reddening - foreground ISM. This study re- 
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Fig. 1. The positions (black and white dots) o f the 89 mainly B-typ e members of Upp Sco are shown (with HD numbers) on top 
of a 100 /urn infrared dust map of this region dSchlegel et al.lfl998l) . The north-east arrows are 3° in length. The pOph cloud can 
be identified by the bright filamentary emission located just left to the center of the figure. The dust emission (on a logarithmic 
grey scale) is proportional to the reddening E^b-v) of sightlines penetrating these clouds; the lowest intensities correspond to E^b-v) 
~0.02 mag (black) and highest intensities to > 2 mag (white). The well-known targets <x Sco (HD 147165, in the pOph cloud) and 
£ Oph (HD 149757, top left) are included in our study. 



vealed small-scale structure - on parsec scales - in the warm neu- 
tral and weakly ionised medium of the Disc-Halo interface. The 
observed low 5797/5780 DIB ratio was found to be consistent 
with the relative high UV radiation levels typically inferred for 
the extra-planar warm medium. 

Nearby OB associations host many bright early-type stars 
confined in a relatively small area of the sky. These stars have 
only few stellar lines in the optical spectrum contaminating the 
interstellar spectrum. Thus, these associations provide a setup 
that is perfectly suited to study the effect of varying local condi- 
tions on the DIB spectrum. One of these associations, Scorpius 
OB2, is a young (5 - 20 Myr), low-density ( » 0.1 M pc~ 3 ) 
grouping of stars divided in th ree subgroups dde Zeeuw et alJ 
Il999t iKouwenhoven et afll2005h . Scorpius (Upp Sco) region is 
the subgroup near the Ophiuchus star fo rming region and the p 
Oph cloud at a distance of 145 ± 2 pc (Ide Zeeuw et al.|[T"999h . 
Combinin g 2MASS extinc t ion m aps with Hipparcos and Tycho 
parallaxes, Lombardi et al.l d2008l) found a distance of 1 19+6 pc 
for the p Ophiuchi cloud (with the core at 128+8 pc~). lMamaiekl 



J20081) suggested a mean distance of 139+6 pc for the distance of 
the Ophiuchus molecular cloud, which they placed within 1 1 pc 
of the centroid of the Upper Scorpius subgroup. 

Filamentary - interstellar - material connected to the 
p Ophiuchus clo ud complex is observed towards Upper Scorpius 
(Ide Geusl Il992l) . The densest part of this complex is the 
pOph dark cloud, a si t e of ongoing low-mass star formation 
dGrasdalen eTal] 119731: iGreene & Yound Il992t IWilking et all 
119971 : iPreibisch & Zinneckerl 1200 ll) that is exposed to the ra- 
diation fields and stellar winds produced by nearby early- 
type stars. A detailed review on the stellar population and 
star formation history of the Sc o OB2 association is given by 
IPreibisch & Mamajekl d2008l) and lWilking et all d2008h . 

The advantages of studying the properties of DIBs in the 
Upper Scorpius region are numerous. It is in close vicinity and it 
has been extensively studied in the past. Detailed information is 
available on both the stellar content (spectral types, photometry, 
distances, kinematics, etc.) and the conditions of the surrounding 
interstellar medium (dust emission and absorption, IR-to-far-UV 
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extinction curves, UV emission, molecular content, etc.)- It ex- 
hibits a significant variation in local environmental conditions 
which should translate into changing properties of the DIBs (if 
they depend on these conditions) when probing different parts of 
the Upp Sco region. 

Previous studies of interstellar gas and dust in the Upp Sco 
region focuse d on the p Oph cloud and a few other nearby 
bright B stars. ISnow et al.l(l2008l) give a concise summary of dif- 
ferent studies of the Upp Sco region covering a range of top- 
ics including UV extinction, atomic and molecular hydrogen, 
atomic and molecular gas, astrochemistry, and DIBs. For ex- 
ample, H2 observations show that sightlines in this region have 
both low (<0. 1) and high (~0. 3-0.6) molecular fractions /h, (e.g. 
ISavage et al.lll977h . 

In this paper we investigate the behaviour of five well-known 
DIBs (at 5780, 5797, 6196, 6379, and 6613 A) and the molecular 
lines of CH, CH + , and CN in the sightlines towards 89 B-type 
stars in the direction of Upp Sco (Fig. [TJ. These targets, within a 
field of 20°x20°, provide a unique and detailed view of the gas 
and dust in this nearby association. In Sect. [2] we introduce our 
sample and provide information on the reduction of the obtained 
spectra. Sect.[3]briefly discusses line-of-sight reddening and dust 
towards Upp Sco. In Sect. |4] we present the observational details 
of atomic and molecular lines as well as diffuse bands. We ex- 
plore the results in Sect. [5] where we discuss first the relation 
between DIB strength, the dust tracer E^b-v), and the molecular 
content. Then, we demonstrate that the DIB ratio 5797/5780 may 
be useful to distinguish between lines-of-sight probing diffuse 
cloud edges and those penetrating denser cloud cores. The ob- 
served differences in physical properties of both types of sight- 
lines are often attributed to the skin effect, the increase in ef- 
fective shielding of molecules f rom UV radiation as one moves 
deeper into an interstellar cloud dKrelowski & Westerlundl 19881 : 
lHerbidfl995l: ICami et all 1 1997b . Furthermore, we have studied 
the effect of local environmental conditions, such as density and 
UV field strength on DIB strengths and ratios. The line strengths 
of CH, CH + , and CN can be used to characterise the physical 
and chemical conditions in the respective sightlines. We have 
constructed a simple dust cloud model to derive the intensity of 
the interstellar radiation field (ISRF) from the observed CH and 
CN line strengths. The paper concludes with a summary of the 
main results (Sect. [6j. 



2. Optical spectra of B-type stars in Upp Sco 

The observed targets cover a region in the sky of approximately 
20°x20° (i.e. -50 x 50 pc at a distance of 145 pc). Within this 
relatively small region 89 sightlines are measured and analysed. 
The positions of the observed targets in the Upper Scorpius re- 
gion are shown on a 100 fim image (far-infrared dust reddening 
map; Fig.[TJ, a reprocessed composite of the COBE/DIRBE and 
IRAS/ISSA maps with the zodiacal foregro und and confirmed 
point sources removed dSchlegel et al.lll998l) . 

Seven out of eight stars that generate 90% of the local inter- 
stellar radiation field (ISRF) are located in this region. The dust 
in these lines-of-sight imposes E(b-v) values from ~0.02 up to 
~0.99 mag (see Sect. [3}, implying local variations in the cloud - 
column or volume - density and structure, and subsequently the 
attenuation of the ISRF. Therefore this large dataset is extremely 
valuable to investigate the effects of environmental conditions 
on the DIB carriers on a local scale. 

Echelle spectra were obtained with the FEROS instrument 
on the ESQ 1.52m telescope at La Silla from 26-30 April 1999. 



T 




.(01) 



Distance [pc] 

Fig. 2. The colour excess E^b-v) (mag) is plotted agai nst the 
distance of the observed targets dPerrvman et al.l H997I) . Stars 
with distance errors larger than 50 pc are shown as squares in 
grey/blue. The low E(b-v) below 120 pc indicates that there is 
little foreground material in front of Upper Scorpius. The in- 
creased scatter at 140 ± 20 pc reflects the density variation as- 
sociated primarily with the p Oph cloud complex. Beyond that 
distance, no substantial increase of reddening is observed up to 
400 pc. A similar figur e for the p Oph cloud region was shown in 
lLombardi et al.l (l2008h (also based on optical photometry, with a 
partial overlap in the selected sky region). 

The spectra were taken at a resolving power of R a; 48, 000 cov- 
ering a spectral range of 3800 to 8500 A. The data were reduced 
using the FEROS context within the ESO-MIDAS data reduc- 
tion package. Data reduction was performed in a standard fash- 
ion, the CCD images were first bias subtracted and subsequently 
the Echelle orders were straightened, extracted, unblazed (flat 
fielded), rebinned (wavelength calibrated) and finally merged. 
We extracted and normalised the spectral ranges of interest. 
Final spectra have signal-to-noise ratios between 100 and 400 
in the wavelength regions of the measured lines. Typically, the 
S/N values are lower for the blue region (-3900-4500 A) with 
respect to the red part of the spectrum. Furthermore, exact values 
differ for each line-of-sight (due to differences in exposure time, 
visual magnitude of the star and the weather conditions). The 
S/N is reflected in the reported equivalent width uncertainties. 

Table Q] (Online) summarises the basic data for the observed 
targets: Henry-Draper (HD) and Hipparcos (HIP) number, spec- 
tral type, right ascension & declination, colour B - V, intrinsic 
colour (B - V)o, reddening E( B -v) (see Sect. |4]», total-to-selective 
visual extinction Ry, and Hipparcos distance (pc). 

3. Line-of-sight reddening and dust 

For each target we derive the reddening E(b-v) from the B and V 
photometry (taken from Tycho-2 and converted to Jo hnson sys- 
tem) after assigning the intrinsic colour (B-V)o (Trom lFitzgeraldl 
[1970) according to the spectral type of the target (as provided 
by the Michigan Spectral C atalog of HD stars; iHoukl 119821; 
iHouk & Srnifh-Moorelll988l) . Visual inspection of the spectral 
range from 4000 to 5000 A, used for th e classification of OB- 
type stars (IWalborn & Fitzpatrick|[T990l) . gives results in good 
agreement with the spectral types listed in the Michigan Spectral 
Catalog. The adopted magnitudes and results for E^b-v) ar e listed 
in Table [T] The total error for E^b-v) is ~0.03 mag, which is 
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Fig. 3. Comparison of the amount of dust, as indicated by E(b-v), 
inferred from optical photometry observations and stellar spec- 
tral classification wit h those inferred from infrared dust emission 
(IRAS 100 fim map: ISchlegel et al.lll998h . p Oph sightlines are 
indicated separately as well as 4 sightlines connected to an in- 
frared dust emission structure north of p Oph. The dashed line 
indicates the 1-to-l relation between the two measurements. 



derived from the error of the Tycho-2 B and V photometry 
(~0.02 mag), the assumed uncertainty (~0.01 mag) in the trans- 
formation to the Johnson system, the colour range of spectral 
sub-types, and the uncertainty in the spectral classification (both 
~0.01 mag for our B stars). 

Two structures of interstellar medium are observed towards 
the Upp Sco complex. A recent study of the distribution and mo- 
tions of the interstellar gas in the p Oph region provides evidence 
for a low density/extinction ISM component, at a distance of 50 - 
80 pc, located in front of the p Oph complex dSnow et al E0081) . 
This nearest sheet-like structure was also observed at a distanc e 
of ~60 pc towards the Sco-Cen region by ICorradi et al.l (120041) . 
This structure has a very low column density and an almost neg- 
ligible effect on the observed reddening. 

The second structure is located at a distance of ~ 1 10 

- 150 pc, consisting of diffuse e xtended portions of the 
dense pOph cloud at 122+8 pc dSnow et all 120081) . This 
is consistent with a mean thickness of ~30 pc found by 
iLombardi et all (120081). Comb i ned w ith the gas densities mea- 
sured by IZsargo & Federmanl d2003h . which suggest a cloud 
thicknesses between 1 and 15 pc, this implies that these clouds 
are not spread homogeneously throughout the Upp Sco region 
but form a patchy complex of scattered and loosely connected 
clouds. 

Note that the column density of this dust sheet (N(H) ~ 3.2 

- 50 xlO 20 cirr 2 : lBohiin et alll978tlDiplas & Savagell 19941) lo- 
cated at approximately 125 pc is an order of magnitude higher 
than that of the nearer sheet. For additional information on in- 
terstellar materia l observed towards the p Oph molecular cloud 
complex see also lMotte et al.l d!998l) . 

The above is supported by the measured colour excess E(p-v) 
as a function of target distance (Fig. |2|. Six stars are probably 
in front of Upp Sco while about 10 to 15 of these are back- 
ground stars. The strong increase of the reddening around a 
distance of 140 pc suggests that most material contributing to 
the extinction is associated with the Upp Sco complex, with 
the observed scatter resulting from variations in the - column 
or volume - density within this region. The re ddening values 
extracted directly from the dust reddening map dSchlegel et alJ 



Il998h are compared with the E(b-v) values obtained from opti- 
cal photometry and stellar classification for the individual sight- 
lines (Fig. [3]). We note, however, that extinction maps based on 
infrared emission or optical/near-infrared star counts show sys- 
tematic offsets with respect to each other and are unreliable at 
small scal es (< 5'), with typi cal lcr uncertainties of 1.2 mag 
in Ay (e.g. ISchnee et al.ll2005T) . For all sightlines the reddening 
inferred from far-infrared emission is higher than that derived 
from optical photometry and spectroscopy, which suggests that, 
for most sightlines, the infrared emission also traces dust that is 
located behind the observed star. Some caution is required com- 
paring these two results as some variation would be expected 
due to calibration/systematic and statistical errors on both val- 
ues. Nevertheless, most sightlines with E(B_y) opt i ca i < 0.3 mag 
have E(B-v)i n frared < 0.6 mag, which is fully consistent with a 
dust sheet associated to the Upp Sco OB association, with both 
stars and dust inter-dispersed with each other. For sightlines in 
the direction of the dense p Oph cloud much higher values for 
E(s-v) (>1 mag) are inferred from the dust map with respect to 
the optical photometric data. Logically, stars visible in this direc- 
tion are likely situated at the front side of this dense cloud (the 
stars at the back will be much fainter / invisible due to higher 
extinction). 

The contribution from foreground material to the observed 
total reddening is very small (E(b-v) ~ 0.02 mag). Therefore, 
we conclude that the dust distribution inferred from both the 
100 jum infrared emission and the line-of-sight reddening is pre- 
dominantly due to the Upper Scorpius complex. In other words, 
the low-density foreground dust sheet contributes very little to 
the total observed values for the infrared emission and the red- 
dening. 



4. Interstellar absorption lines 

the interstellar absorption lines observed towards the 89 Upper 
Scorpius targets. We determined equivalent widths for the five 
strong DIBs at 5780, 5797, 6196, 6379, and 6613 A, for the di- 
atomic lines of CH, CH + , and CN as well as for the Ki and Cai 
lines (Table|2|. The Nai D doublet is omitted because it is satu- 
rated for the majority of sightlines. Line profiles and central he- 
liocentric velocities for the atomic and di-atomic lines are given 
in Fig. ID. II and Table [5] (Online), respectively. To illustrate, the 
velocity absorption profiles of the atomic and molecular absorp- 
tion lines towards HD 147889 are shown in Fig. [4] 



4. 1 . Diffuse interstellar bands 

Although more than 300 DIBs are known we focus here on the 
five strong and narrow bands at 5780, 5797, 6196, 6379, and 
6613 A. The strength and width of these features facilitates the 
measurement of modest column densities of their carriers in 
slightly reddened sightlines. Additionally, the Galactic relation- 
ships between DIB strength and reddening are well established 
for these DIBs both in the Galaxy and beyond. The equivalent 
width, W, is measured via a straight line continuum integration 
across the absorption feature (see Appendix. lAl Online). For the 
DIB measurement we do not expect significant contamination 
from stellar atmosphere lines (see Appendix. |Bj Online). The 
measured equivalent widths, or 2<x upper limits, are listed in 
Table |2] (Online) for the five DIBs towards the 89 targets in Upp 
Sco. This is the first consistently measured data set containing 
this many sightlines within one region. 
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Fig. 4. Example of the velocity profiles of interstellar absorption 
lines. The line-of-sight shown is that towards HD 147889. From 
top to bottom the CH + (both at 3957.7 and 4232.5 A, latter as 
thick line), CH, CN R(0), Ca 1, and K 1 line (both 7665 and 7699, 
latter as thick solid line). The strong feature at ~15 km s in 
the bottom panel is a telluric absorption line. Similarly, velocity 
profiles for all sightlines are included in Appendix iDl (Online) . 
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To ascertain the accuracy and consistency of our results 
we compa re our measure d values with those a v ailable in the 
litera tur e (iHerbid 119931: " 



200 U IGalazutdinov etal] l2002t 



Seab & Snowl 119951: 



Thorburn et 



Megier etail 
alJ 12003b 



Galazutdinov et al.l I2004L iMegier et all l2005t ISollerman et alJ 
2005t and Friedman et al.ll201 lh . The result is shown in Fig. [5] 
for the 5780 and 5797 DIBs, and in Fig. (Online only) for the 
6196, 6379, and 6613 A DIBs. The correlation between values 
from this work and from literature are good (for the 5 DIBs, 
the correlation coefficients r range from 0.92 to 0.99). The 
linear regressions reveal a small non-zero offset, indicating that 
our values are systematically lower by a few percent. We note 
that for individual cases the values between different studies 
vary significantly by as much as 20%. Small inconsistencies 
(between all studies) arise naturally from differences in the data 
quality (S/N, resolution) as well as differences in the adopted 
methods for equivalent width measurements (adopted stellar 
continuum, contamination from nearby weak features, adopted 
integration limits, inclusion or removal of underlying broad 
bands). In conclusion, the measured equivalent widths are 
consistent with previous studies but do show a small systematic 

offset. 

ICami et al.l (fl997h found that the AA5797, 6379, and 6613 
DIBs show a good correspondence to each other, with r ~ 0.8. 
These authors also found that the /15780 DIB is moderately cor- 
related with the ,16613 DIB (r - 0.65) and weakly to the 6379 
DIB (r = 0.47). Recently, iFriedman etaP d201 if) found high 



Fig. 5. Equivalent widths of the /U5780 and 5797 DIBs 
previously measured for sightlines included in this survey 
compared to values measured in this work. For several 
sightlines several literature values are available, illustrat- 
ing the "intrinsic" scatter in equivalent widths due to 
measur e ment m e thods. Literat ur e val ues are taken from 
IHerbid dl993l): ISeab & Snowl (Il995h: IGalazutdinov et all 
d2003l): IGalazutdinov et all d2004l): 



d2002b: iThorburn et al 



Sollerman et all d2005h : iMegier et all d200ll l2005h : and 



iFriedman et alJ (1201 lh . Linear regressions are shown in each 
panel. Correlation coefficients r are 0.99 and 0.96 for the 5780 
and 5797 A DIBs, respectively. Slopes and intercepts of these 
regressions are 1.02 and 11 mA for the 5780 A DIB and 1.12 
and -6 mA for the 5797 A DIB, respectively. 



values for r (ran ging from 0.93 to .99) for the Galactic DIB 
pairs in Table [3] McCall et alJ d2010h reported a nearly perfect 
correlation (r = 0.99) between the ,U6196 and 6613 DIBs to- 
ward 1 14 Galactic diffuse cloud sightlines. In this work we find 
r = 0.8 for /LI6196-6613 DIB pair, which is less than for other 
pairs. The correlation coefficient between the five DIBs mea- 
sured in this study are given in Table [3] In line with previous 
results, the ,15797 DIB has a good correlation with both /L16379 
and 6613 DIBs, however, it shows a poor correlation with both 
AA5780 and 6196 DIBs. The ,1,16379 and 6613 DIB pair shows 
the strongest correlation, with r = 0.92. In fact, the ,16613 DIB 
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correlates well with all four DIBs. The /15780 DIB shows a good 
correlation with the /16613 DIB (r = 0.85). Unexpectedly, the 
other two DIB family members of i6613 (i.e. AA5797 and 6379 
DIBs) have a weaker correlation. Re stricting the computatio n of 
r to the 13 sightlines present in both lFriedman et all (1201 ll) and 
this work, increases r for our data (but lowers r slightly for the 
Friedman sample). For example, for the /LI5780-5797 DIB pair 
r = 0.75 (this work) and r = 0.86 (Friedman); for the /L16196- 
6613 DIB pair r = 0.97 (this work) and r = 0.99 (Friedman), and 
for the ,M5780-6196 DIB pair both studies give r = 0.93 (for the 
complete dataset Friedman report r = 0.97). The higher Pearson 
correlation coefficient (independent of quoted error bars) for the 
Friedman data suggests that the overall uncertainties on the mea- 
surements are lower than for this work, resulting in an improved 
correlation. Partly this is due to the fact that our sample includes 
a large fraction of sightlines with low values for E(b-v), and thus 
weak DIBs. On the other hand, restricting the comparison to the 
Upp Sco sightlines in common lowers r in both samples (prob- 
ably as there are fewer data points), and also reduces the dif- 
ference between the two sets. This could be partly due to an 
increased effect of local variations in the DIB spectrum on the 
correlation coefficient (such effects would be averaged out in 
a larger Galactic survey probing many different regions as op- 
posed to probing a peculiar region like Upp Sco). 



Table 3. Pearson correlation coefficients r between the observed 
DIBs. 



DIB 5780 


5797 


6196 


6379 


6613 


5780 1 


0.72 


0.74 


0.75 


0.85 


5797 


1 


0.69 


0.87 


0.85 


6196 




1 


0.81 


0.80 


6379 






1 


0.92 


6613 








1 



4.2. Molecular lines 

We have measured equivalent widths and heliocentric radial 
velocities for the CH (^- est = 4300.313 A), CH + (A If , st = 
4232.548 A), and CN R(0) (A KSt = 3874.608 A) lines (Tables 
and [5] respectively; Online). In a few (about 5) cases, the CN 
(3874.608 A), CH (3886.410 A), and CH + (3957.70 A) lines are 
tentatively detected (see e.g. Fig. These lines are weak and 
have large (>50%) uncertainties. It may be that the strongest 
CH line is saturated, which can occur for individual compo- 
nents with W(CU) > 20 mA (IVan Dishoeck & Blacklll989l) . The 
CN R(0) transition is also prone to saturation for (individual) 
components with W > 6 mA, leading to underestimated column 
densities, though the corrections are less than about 20% up to 
W = 15 mA dSlvk et al.l2008l) . The CN lines toward HD 147683, 
HD 147701, HD 147889, and HD 147932 likely suffer from sat- 
uration. For the lines-of-sight including the strongest CN lines 
in our sample, HD 147932, HD 147701, and HD 147889, the 
column densities would need to be corrected by a f actor 1.27, 
1.7, and 2.1, respectively (following Slvk et al. 200jl adopting a 
value of 1 km s _1 for the Doppler broadening). Saturation also 
occurs for CH + if W >20 o r >40 mA for components with b — 1 
or 2 km s _1 , respectively (I Allenl IT994h . Only a few sightlines 
have measured total W larger than these limits, and even for 
these cases the individual (unresolved) velocity components are 



not expected to be strongly saturated as noted above. For the 
sightlines towards HD 147683, HD 147889, HD 147933, and 
HD 149757 the equivalent width ratio between the (tentatively) 
detected weaker and stronger lines of both CH and CH + are close 
to - within the uncertainties - the expected ratio of ~3.9, and 
~1.9, respectively. Only for the latter two sightlines are these 
ratios significantly lower (~1.5) indicative of some saturation. 

For certain spectral types stellar line contamination can 
complicate measurements. However, for the majority of spec- 
tra presented here this problem could be well resolved (see 
also Appendix |Bl). Equ i valent widths for CH and CH + given by 
iFederman et alTTl994h : lMegier et al.l (l2005l) . and lWeselak et al.l 
d2008bl) for nine Upp Sco sightlines in common with this work 
are consistent with our values. Furt hermore, the repor ted W(CN) 
are consistent with values given bv lStvk etaU (120081) for six tar- 
gets in common, although their reported error bars are smaller. 
The velocity profiles for CH, CH + , and CN are also shown, for 
the relevant sightlines, in Figure [Dj] (Online). 

4.3. Atomic lines 



Inspection of the Nai (5889.951 & 5895 924 A ; lMortonll2003l) 
and Ki (7664.91 & 7698.974 A; lMortonll2003h doublets shows 
that most (75 of 89) sightlines are dominated by one strong ve- 
locity component. This strong component displays asymmetries 
and broadening for a number of sightlines suggesting that in real- 
ity multiple unres olved narrow com ponents may be present (see 
e.g. Sect. 15.51 and ISnow et al.ll2~008h . The obtained spectral res- 
olution is not sufficient to resolve hyperfine splitting (order of 
~1 km s" 1 ). For 10 sightlines (all with low reddening; E(b-v) 
< 0.2) two or three weaker components, clearly separated in ve- 
locity space, could be discerned in the Nai profiles. For a small 
number of sightlines we also detect Ca i at 4226.73 A. Equivalent 
widths, heliocentric radial velocities, and profiles for Ki and 
Cai are included in Tables|2]and|5] and Figure ID~T1 respectively 
(Online). Equivalent widths and profiles are not provided for the 
highly saturated Na i doublet as these preclude any column den- 
sity measurements. However, approximate central velocities are 
included in Table[5](Online). Similarly, the Ca n line is saturated, 
but also suffers from stellar contamination and reduced spectral 
quality in the blue. 

5. Results & Discussion 

In the following we present and discuss the relation between 
equivalent widths of the observed DIBs, molecules, and the 
line-of-sight reddening. We investigate whether the correlations 
we found can be explained in terms of the skin effect and ex- 
plore the spatial variation of DIB strength and strength ratios. 
Furthermore, we discuss the velocity structure of the ISM, as 
well as a model of the dust sheet and the inferred effective inter- 
stellar radiation field (ISRF). The ISRF strength, /uv, and molec- 
ular hydrogen fraction, /h 2 , are both discussed in view of cr and 
(-type clouds. 

5.1. DIBs and dust 

In Fig. we show W^(5780) (top) and W A (5197) (bottom) 
against E^-vy Both the average Galactic and Upp Sco relation- 
ships are shown. Several conclusions can be drawn immediately 
from this initial result: 1) A linear model does not adequately 
describe the relation between the measured values (x 1 » 1); 
2) The average DIB strength per unit reddening in Upp Sco is 
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similar to the Galactic average; 3) There is a positive trend be- 
tween the amount of DIB carriers and the amount of dust in 
the diffuse ISM; 4) There is a significant scatter from this mean 
linear relationship (which is also observed for the Galaxy-wide 
surveys), especially for the /15780 DIB. This results in a poor 
Xred- I n particular, for sightlines with E ( b-v)~ 0.2 to 0.3 mag 
(which would typically be expected to be single diffuse cloud 
sightlines) there is marked range in strength of the DIBs (for 
both /L/L5780 and 5797 DIBs the strength can vary by factor of 
about four to five). The scatter (standard deviation) around the 
mean is equally high for higher E(b-v), but for those multiple 
cloud components are more likely to contribute and confuse the 
true variations in individual clouds. The strength-reddening rela- 
tions for the AA6196, 6379, and 6613 DIBs are similar to that for 
/15797, albeit with different slopes and an increased scatter (see 
Fig- [SJ- The linear fit method using uncertainties in both param- 
eters is a n implementation of the routine fitexy from Numerical 
Recipes dPress et al.ll 19921) where = x 2 /(N - 2), with N the 
number of data points. A good fit will have^-^ « 1. Despite 
significant intrinsic scatter in the DIB versus reddening relations 
the least-square linear fit results are given in Table |4]to facilitate 
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Fig. 7. Equivalent width versus 
5797 A (bottom) DIBs . The average 

(dashed: ICox^td]|200l and dotted: Friedman et alj|20lTb and 150 T X rert = 6.4 

the linear least-squares fit for Upp Sco (solid; Sect. [5} are shown. 
Reduced chi-squared values 0^ ed ) for the latter are indicated in ^ 
the respective panels. Intercepts and slopes for the linear least- g 
squares fits are given in Table|4] To avoid biases, the regressions ^ 
were not forced to go through the origin, and upper limits for the ^ 
Upp Sco data were not taken into account. S- 
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Fig. 8. Equivalent width versus E ( b-v) for the 6196, 6379, and 
6613 A DIBs observed towards the Upp Sco lines-of-sight. The 
reduced chi-squared OkfjL/) f° r the linear fits (dashed lines) are 
indicated in the respective panels. The linear fit parameters are 
given in Table [4] 



estimates of interstellar line-of-sight reddening from observed 
band strengths. 

These deviations could reveal the effects of local condi- 
tions on the balance between DIB carrier formation and de- 
struction (including changes in e.g. ionisation and hydrogenation 
state), and therefore the abundance and physical properties of 
the DIB carrier. The generally positive correlation between DIB 
carriers and reddening suggests a link between the presence of 
dust grains and the molecules responsible for the diffuse bands. 



8 



Vos, Cox, Kaper, Spaans & Ehrenfreund: Diffuse interstellar bands in Upper Scorpius 



Table 4. Slopes and intercepts of the linear least-square fits. 
Uncertainties in both coordinates are taken into account. Non- 
detections and upper limits were excluded from the fit proce- 
dure. The fits were not forced to go through the origin, though 
it can be noted that in most cases the derived intercept is within 
2<x of the origin. These relations can be used to derive estimates 
for the (interstellar) line-of-sight reddening from measurements 
of the diffuse band strengths. 



Correlated parameters 


Intercept 


Slope 


X-red 


r 




(mA) 


(mA/E (2! -v)) 






W(5780) - E B _ V 


2.7 ± 3.3 


462.0 ± 12.7 


10.8 


0.79 


VK(5797) - E B - V 


-5.0 ± 1.1 


159.0 ±4.1 


4.0 


0.92 


W(6196)-E B _ V 


3.8 ± 0.4 


35.2 ± 1.8 


6.0 


0.72 


W(6379) - E B - V 


-1.5 ±0.8 


88.5 ± 0.8 


5.6 


0.85 


W(6613)-E B _ V 


-2.2 ± 1.2 


177.8 ± 4.6 


6.4 


0.86 


^(5780)^-^ 


3.0 + 7.8 


640.2 ± 43.6 


3.3 




W(5780) f - E B -v 


-23.8 ± 5.0 


419.2 ± 13.7 


9.3 




W(5797) lr -£ B _„ 


-1.7 ± 1.8 


127.0 ± 9.6 


2.8 




W(5797)<- - E B - V 


-0.5 ± 1.9 


153.3 ±5.3 


5.1 





Figs. IC. 1 l and lC.2l (Online appendix) illustrate that there are large 
variations, particularly at intermediate E(b-v)~ 0.2-0.3 mag, in 
the DIB strengths normalised by the amount of dust in the sight- 
line. At lower E^b-v) the measurements are inaccurate, and at 
higher E^b-v) the presence of multiple clouds in the line-of-sight 
appears to reduce the effect of variations in individual clouds on 
the composite, total line-of-sight DIB spectrum. The behaviour 
of the DIBs in relation to molecular tracers and the local envi- 
ronmental conditions will be discussed in the next sections. The 
different behaviour of the /M5780 and 5797 DIBs is used as a 
tool to study the deviations of both DIBs from the mean trend 
with E (B _v). 



5.2. The skin effect 

DIB carriers seem to reflect the evolutionary cycle of molec- 
ular carbon species (such as aromatic molecules) t hrough for- 
mation, ionisation, recombin ation, and destruction dCami et al.l 
ll997tlRuiterkamp et al]l2005l) . 

Uncharged aromatic molecules exhibit strong absorption 
bands in the UV and visible (blue) range while their cations and 
anions show s pecific transitions in the visible (green-yellow) and 
near-infrared ( ISalama et al.lll999l) . Each DIB carrier is thus in- 
fluenced by the interstellar radiation field in a particular way, 
since its molecular properties such as ionisation potential and 
electron affinity are unique. 

Interstellar clouds are exposed to t he interstellar radiatio n 
field which drives their photochemistry (ISnow & McCallll2006r) . 
The UV radiation is attenuated (by dust) increasingly from cloud 
edge to core, giving different steady-state solutions for the pho- 
tochemical reactions (like the ionisation-state) in different parts 
of the diffuse cloud. Thus interstellar species are subjected to 
stronger radiation at the edge than in the centre of the cloud. 

Especially the (molecular) DIB carriers are believed to be 
sensitive to UV radiation. The signatures of more stable DIB car- 
riers (such as corresponding to the /15780 DIB) show a relative 
higher intensity in lower density, higher Tuv regions with respect 
to less stable DIB carriers which are more rapidly destroyed at 
high 7uv {e.g. A5191 DIB). These reach higher intensity only 
in more UV protected denser regions (where more stable DIB 
carriers like the /15780 DIB are less efficiently ionized and thus 
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Fig. 9. W,,(5780) (top) and W,,(5797) (bottom) versus E (B -v) for 
the cr and ( subgroups, respectively. Error bars, identical to those 
in Fig. [7] have been omitted for clarity but taken into account for 
the linear fit. Parameters for the linear fit and regression are given 
in Table |U 



reduced in strength). For a sightline probing Upper Scorpius a 
larger amount of dust is expected to correlate with, on average, 
higher densities, especially as there often is only one apparent 
strong interstellar velocity component. 

The effect of shielding (to a certain degree) of molecules 
from strong UV radiation is often referred to as th e skin-effect 
(e.g. iKrelowski & Westerlundlll988t lHerbiglll995h . The skin- 
effect reflects the life cycle and charge distribution of DIB carri- 
ers, which can can lead to an interpretation of high DIB carrier 
concentrations in the outer cloud layers. However, DIB carriers 
are also expected to be present in high concentrations in denser 
regions although in a different charge state (neutral) that can only 
be observed in the UV. 

ICami et all (1 19971) inferred that the A5180 DIB carrier 
reaches its maximum abundance when exposed to the interstellar 
UV radiation field (typically near the edge of a cloud), whereas 
the A5191 DIB carrier is more easily ionised and destroyed. Even 
more, at very low E(b~v) (<0.1 mag) only ver y few DIB carri- 
ers su rvive due to the high rate of UV photons ( Jenniskens et alJ 
119941) . The relative abundance between the AA57S0 and 5797 
DIBs reflects an interplay between neutral, ionised, and de- 
stroyed DIB carriers along the entire line of sight. This balance 
is affected not only by the impinging radiation field, but also 
by the carbon abundance and the dust particle size distribution 
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Fig. 10. The W(5797)/W(5780) ratio plotted against E (B -v)- The 
distribution peaks at an E(b-v> of ~0.25 mag, indicating most op- 
timal conditions for formation of the 5797 A DIB carrier and the 
destruction c.q. insufficient excitation of the molecule giving rise 
to the 5780 A DIB. Nonetheless, the significant scatter suggests 
that processes additional to dust extinction are important. 



(ICox & Spaansl2006h . A difference in the observed ratio of these 

two DIBs is th u s direc tl y related to the skin-e ffect. 

Krelowskil dl989t) : ISneden et alJ dl99lh : IKrelowski etafl 



(fl992¥ identified two types of clouds, referred to as cr and f- 
type, cr-type clouds show atomic lines and DIBs, but the molec- 
ular lines are weak or absent, while f-type clouds have strong 
diatomic lines in addition to DIBs. The main difference be- 
tween both types lies in a combination of density and UV ir- 
radiation by the ISRF, with cr clouds associated with low den- 
sity and/or strong exposure to UV radiation, while f clouds are 
associated with higher densities and/or more protection from 
UV radiation. Therefore, differentiation between cr and £-type 
clouds is directly linked to the skin-effect described previously. 
For sightlines probing ^-type clouds the 5797 A DIB is deeper 
than the 5780 A, while for cr-type clouds the reverse is ob- 
served. Therefore, the W(5797)/W(5780) ratio has been used 
to distinguish between UV exposed (cr) and UV protected {() 
sightlines. The nomenclature for the cr and f type sightlines 
is historical and based on the repres entative lines-of-sight to- 
wards cr Sco and C Oph; respectively (IKrelowski & Westerlundl 
119881; IKrelowski et al'l 1 19921 IKrelowski & Snedenlll995l) . Note 



that both sightlines are included in our analysis. 

In this work we re-establish this classification, assuming a 
relatively equal distribution of sightlines probing dense versus 
diffuse clouds. Sightlines are classified cr when the ratio is lower 
than the weighted mean of the ratio minus lcr, while ratios 
higher than the weighted mean plus lcr are classified as The 
remaining sightlines are classified as intermediate. The results of 
this selection for individual lines-of-sight are included in Table|2] 
(Online). Figure [9] shows that the application of our classifica- 
tion to the data in Fig. Qimproves the relation (reduced scatter) 
between DIB strength and reddening. Indeed, Fig. [9] shows that 
W(5780)-E(b_i/) has an improved reduced x 1 f° r the cr and (, 
sightlines respectively, though only a marginal improvement is 
found for W(5797)-E (B „ V) (where higher W(5797) tend to corre- 
spond to £-type environments) revealing the /15780 DIB is pri- 
marily giving rise to variations in the W(5797)/W(5780) ratio. 

Note that the original classification is based on the central 
depth, A, of the two DIBs. If A5797 > A5780 (corresponding to 



W(5797)/W(5780) > 0.4) the line-of-sight is considered as 
type. Increasing our selection threshold for the DIB ratio to 0.4 
would imply that the "intermediate" sources would be included 
in the cr-group as well as a few ^-types. However, sightlines with 
both low (/h,<0.3) and high (/hi>0.4) molecular content cur- 
rently classified f would also be re-assigned as cr-type. We note 
that it is impossible to make a sharp distinction between cr and 
f-type sightlines as there is - as expected - a smooth transition 
of physical conditions characterising both types. 

Fig. [10] shows W(5797)/W(5780) as a function of redden- 
ing. The distribution peaks at an E^b-v) of ~0.25 mag, indicating 
most optimal conditions for formation of the 5797 A DIB carrier 
(sufficient shielding), or alternatively less optimal conditions for 
the carrier associated to the 5780 A DIB carrier (insufficient UV 
photons to transform it into its ionic form and thus not absorbing 
at the visible wavelength). For sightlines with E^b-v)>0.4- mag 
the conditions for formation of the /15797 DIB are sub-optimal, 
but still more favourable with respect to the /I5780 DIB car- 
rier than for cr sightlines. At very low E(b-v) (<0.1 mag) the 
/15797 DIB carrier is under-abundant (due to more efficient de- 
struction of molecules by the stronger ISRF) with respect to the 
/15780 DIB carrier. Also, the W(5797)/W(5780) ratio itself dis- 
plays a bimodal distribution with a strong peak at about 0.2±0.05 
(Fig. |24l Sect. I5.71 i. The sightlines associated to this peak are 
predominantly cr-type (which indeed we may consider to rep- 
resent typical diffuse ISM though this should be confirmed by 
studies of other regions). There is a second smaller peak "bump" 
at ~0.45, corresponding to the £-type sightlines. Although, like 
the /15797 DIB the AA 6196, 6379 and 6613 DIBs also show 
signficant scatter on the respective W-E(b-v) trends, the /15780 
DIB reveals the clearest distinction in behaviour between cr and 
£ sightlines. The link between the W(5797)/W(5780) ratio and 
the strength of the ISRF is discussed in more detail in Sect. 15.61 

For comparison we plot also the W(6196)/W(6613) ratio as 
a function of E(b-v) in Fig. IC.4I This ratio is less sensitive to red- 
dening and therefore is not such a useful tracer of local condi- 
tions such as density and UV irradiation. This is indeed expected 
from the recent results bv lMcCall et al.l (1201 Ol) who found an ex- 
cellent correlation between the 6196 and 6613 A DIB strengths. 

5.3. Diatomic molecules and dust 

Different interstell ar species are r estricted to different regions 
(see e.g. Fig. 6 in IPan et"aT1l2005l) : CN and CO are present in 
dense regions, CH and K 1 are predominantly present in moder- 
ately high density regions (« > 30 cm 3 ), and CH + and Ca 1 in 
intermediate density regions (« ~ 10- 300 cm- 3 ). 

In Fig. [IT] we show W(CH) (top) and W(CU + ) (bottom) 
versus E(b-v)- It can be seen that CH correlates much better 
with E( B-y) than CH + , which is in line with previous obser- 
vations dCrawfordlll989l IKrelowski etaD 1 19991) . CH traces the 
dense, molecular gas and its abundanc e is directly prop ortional 
toN(H 7 ) asN(CH)/N(H ? ) = 3.5 x 10 8 (iFedermanl 19821: iMattilal 
119861: IWeselak et all 12004 IShefferet al.ll2008l) . For 8 lines-of- 
si ght direct measureme nts of N(H2) (IUE or FUSE; compiled 
in lFriedman et al.ll201 lb can be compared to those derived from 
N(CH) in this work (Table 0. The scatter i s less than ~0 . 5 dex , 
and in good agreement with the results of She ffer et al.1 (120081) 
and references therein. Theoretically, one can thus infer the 
molecular hydrogen fraction fu 2 from N(H2) derived from CH 
together with N(Hi) d erived from WT5780) flog N(Hi) = 19.00 
+ 0.94 log(M5780)): lFriedman etafll201 1|) . The resulting val- 
ues for /hi are given in Table|2] These values are consistent with 
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Fig. 11. Equivalent width versus E {B - V) for the CH 4300 A (top) 
and CH + 4232 A (bottom) transitions. For <x (black squares) and 
£-type (blue circles) sightlines designations see Sect. 15.21 and 
Sect. GO 



- tho ugh systematically hig her than - the directly measured /h, 
(e.g. iFriedman et alJuOl It Table |2]i for the eight sightlines in 
common. Here we have used the average Galactic relation be- 
tween W(5780) and N(Hi), whereas this relation may actually 
be lower fo r Upp Sco (similar to the lower gas-to-dust ratio in 
this region: ISchlegel et al.lll998h thus leading to a higher esti- 
mate of /h,. The strongly improved regression coefficient be- 
tween CH and E(b-v) for the £-type sightlines (r = 0.83) com- 
pared to its cr-type equivalent (r = 0.53) supports the interpre- 
tation that £-type lines-of-sight trace dense gas. It is noteworthy 
to recall that the significant scatter for the diffuse band strengths 
at low E(£_v) (~0.2 mag) as illustrated in Figs. [7] and [8] is not 
obs erved for CH. CH + , o n the other hand, is not a good tracer of 
H2 dWeselak et al.ll2008ah . Therefore, the low value of the corre- 
lation coefficient for CH + is not unexpected. Furthermore, note 
that si gnificant amounts of C H are needed before CN is pro- 
duced ^FedermanetaL 1984j) 2 _whh_the latter tracing also rela- 
tively dense material (Uoseph et al.lll986l) . 

Work bv lCrawfordl ( ll989l) suggests that the ratio of N(CH + ) 
and N(CH) is indicative of the turbulent or quiescent nature of 
the interstellar medium in the line-of-sight. For shocked envi- 
ronments an offset velocity between CH and CH + or a velocity 
broadening of CH + is predicted by models. For the sightlines 
in this work we obtain an average Av of 0.3 km s , with in- 
dividual velocity measurements that have errors of about 1 to 
2 km s (see also Sect. 15.51 ). Our data support recent surveys 



Fig. 12. W(CH + ) versus W(CH). Assuming optically thin lines 
the N(CH + )/N(CH) ratio is equal to 0.95x W(CH + )/W(CH). 
Turbulent and quiescent regions are indicated by shaded areas. 
See text for further details. For cr and £-type sightlines designa- 
tions see Sect.l5.2land Sect. 15.71 



whic h find no evidence for a velocity diff erence between CH and 
CH + (ICrane et al.lll99l iPan et al.ll2005l) . The data do not allow 
for an accurate measurement and comparison of CH and CH + 
line widths. The line profiles of atomic and di-atomic species 
can be compared in Fig.|4](and associated Fig. ID.ll in the Online 
material). W(CH + ) is plotted against W(CH) in Fig.[12]with the 
turbulent (N(CH + )/N(CH) > 2) and quiescent (N(CH + )/N(CH) 
< 0.5) regions indicated by the shaded areas. T he general cor- 
relation between CH + and CH dPan et al.ll2005l) is poor, but it 
appears that two separate trends might in fact exist for the qui- 
escent and turbulent regions, respectively, potentially indicative 
of different dominant CH + production mechanisms. The dense 
cloud tracer CN is only detected towards £-type lines-of-sight, 
supporting the interpretation that the latter probe dense clouds. 
The cr and £ type sightlines show different trends for W(CH), 
but not so clearly for W(CH + ). 

5.4. DIBs and small molecules 

In this section we discuss the behavi our of the AA5797 and 
5780 DIBs with respect to CH and CH + . lKreIowski et alJdl992h 
found that CH and CN are only detected if the A 5191 DIB is 
deeper than the ,15780 DIB. IWeselak et all d2008bl) studied CH, 
CH + , and CN in relation to DIBs for a large, inhomogeneous 
sample of sightlines. These authors found a good correlation 
between W(5797) and N(CH), but a poor correlation between 
W(CH)/E (B _ y) or W(CN)/E (B -v) versus W(5797)/W(5780). The 
correlation of W(5797) vs. N(CH) is further improved by exclud- 
ing sightlines with overabundant CH. Their conclusion is that 
the /15797 DIB carrier is favoured in environments with higher 
molecular gas content. On the other hand CN traces a denser 
medium where the production of the A5197 DIB is apparently 
more inefficient. 

Figures [T3] and [14] show the relationship between molecu- 
lar line strengths (VK(CH) and W(CU + )) and diffuse interstellar 
band stre ngths (W^(578 0) an d Wj(5797)). These re sults are in 
line with iHerbid (1 19931) and IWeselak etail d2008bl) . who con- 
cluded that DIB strengths correlate better with E(b-v) and Hi 
than with any other feature originating from the gas phase. These 
DIBs have a st ronger correlation wit h CH than with CH + (this 
work) or CN dWeselak et al.l l2008bl) . The positive correlation 
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Fig. 13. W(CH) (top) and W(CH + ) (bottom) versus W^(5780). 
The i5780 DIB shows no direct correlation with either CH or 
CH + . Looking separately at the cr and ^"-type sightlines one can 
distinguish different behaviour between the molecular lines and 
the /I5780 DIB for both types. Sightlines classified as intermedi- 
ate are indicated by green crosses. DIBs are stronger with re- 
spect to CH and CH + line strengths for cr-type sightlines. In 
other words, the 5780 DIB carrier abundance is lower for £-type 
clouds which have a higher molecular content. 



with E(b-v) suggests that even though grains do not give rise 
to the diffuse bands they do play an important role in the ei- 
ther the DIB carrier formation - via e.g. grain surface reactions - 
or destruction - e.g. attenuation of UV radiation - processes. The 
CH molecule and the A5191 DIB correlate tightly, indicating that 
the A5191 DIB carrier is most abundant in CH / H2 clouds (see 
also IWeselak et all 120041) . Some correlation is expected since 
both species correlate with E^b-v)- For individual clouds a larger 
W(CH) is indicative for the formation in denser clouds, which 
explains the tighter correlation with the /15797 DIB compared to 
the /15780 DIB. Note however, that the strongest molecular fea- 
tures potentially arise from (unresolved) multiple components 
of the ISM which are not necessarily denser (see e.g. velocity 
profiles for K 1 in Fig. ID. 11 1. 

Figure [TJ shows the Wi(5797)/Wi(5780) ratio versus the 
W(CH) (top) a nd W(CU + ) (bot t om) n ormalised to E {B - V) . In 
agreement with iKrelowski et al.l d!999l) . a stronger correspon- 
dence is observed for the W / j(5797)/W,i(5780) ratio versus 
W{CH)/E (B -v) compared to that for W,i(5797)/W / i(5780) versus 
W(CH + )/E(B_v). This confirms that ^-type clouds (dense, A5191 
DIB favoured) are connected to a higher molecular content, im- 
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Fig. 14. W(CH) (top) and W(CH + ) (bottom) versus ^(5797). 
There is a moderately good correlation (r = 0.84) between 
W(CH) and W A (5797). There is no correlation (r = 0.47) be- 
tween CH + and W.^5797). 



plying furthermore that the DIB ratio is related to the abundance 
of cold cloud molecular species and properties of interstellar 
dust as suggested in Sect. 15.21 The poorer correlation between 
this DIB ratio and W(CVL + )fE (B -v) then suggests that CH + forms 
in regions with different conditions, such as in the clouds outer 
edge, where the UV radiation field is much stronger. 

In Fig. [16] the W /i (5797)/W /l (5780) ratio is plotted against 
the W(CH + )/W(CH) ratio. Again, turbulent and quiescent ISM 
are indicated. This plot reveals no marked correlation between 
these ratios. Tentatively, it shows a high DIB ratio (i.e. £-type) 
for quiescent clouds and a low DIB ratio {i.e. cr-type) for turbu- 
lent clouds, which supports the idea that both CH and the A5191 
DIB trace moderately dense regions, while CH + traces the cloud 
edges and inter-cloud regions. Although both the /15780 DIB 
and CH + are related to the outer edges of diffuse clouds they 
do not reveal a strong correlation (although there appears to be 
a positive trend when considering only cr-type sightlines) and 
thus appea r to react to ch a nges i n the ISRF differently. In agree- 
ment with IWeselak et all (l2004l) and Sect. 15.21 it seems that the 
Upp Sco region is somewhat turbulent, but is absent of extreme 
shocks. However, uncertainties in W(CH) and W(CH + ) are too 
large to draw firm conclusions. 
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Fig. 17. The heliocentric radial peak velocity for K i and Na i are shown in a sky coordinate (right ascension and declination) plot. 
The radial velocity of the gas (in the diffuse ISM) in Upp Sco is highest in the upper-left corner, the material is approaching with 
velocities up to 15 km s~' (along the line-of-sight). The gas in the lower-right corner has the lowest velocity (with respect to the 
Sun). (See on-line electronic version for colour figure). 



5.5. The ISM velocity distribution 

We measured the heliocentric radial velocities for both atomic 
and molecular lines towards Upp Sco. The strongest interstellar 
lines are observed at a radial velocity of about -9 km s , and a 
weaker absorption component is detected at about -22 km s~' 
(e. g. Fig. [Hi. This is fully in-line with recent results reported 
bv lSnow et ail (120081) who studied 16 lines-of-sight towards the 
Upp Sco region. The velocity component of -9 km s~' corre- 
sponds to the patchy dust sheet at a distance of 110 - 150 pc, 
which is loosely connected to the p Oph dense/molecular cloud 
at 122 pc (see also Sect.0. The other, weaker velocity compo- 
nent (v = -22 km s~') is linked to the tenuous low density dust 
layer at ~50 pc. The average velocity difference between CH and 
CH + is 0.3 km s~' (for the Ki doublet lines the average velocity 
difference is 0.07 km s~'). Within the limits of the observations 
we confirm that there is no evidence for a CH-CH + velocity off- 
set in Upp Sco which is predicted by models for regions with 
strong shocks. 

The relatively broad DIBs preclude a detailed radial veloc- 
ity determination (for the obtained S/N, resolving power, and 
spectral quality). However, first order estimates (for the sight- 
lines with strong DIBs) show no systematic differences between 
atomic, molecular, and DIB velocities. For the relatively narrow 
5797 A DIB we measure radial velocities roughly between -20 
to km s _1 . 

This large set of radial velocity information allows us to map 
the velocity of the diffuse ISM c louds in front of th e observed 
stars, very similar to the work bv lSnow et al.l (120081) . In Fig. [17] 
we show the color-coded interstellar radial velocity of K i and 
Na i as a function of declination and right ascension. Assuming 
that indeed the observed sightlines probe different parts of a sin- 
gle dust-sheet, it thus appears that this sheet is moving differ- 
entially. The upper-left corner of the sheet (where most young 
OB stars are) is moving towards us while the lower-right corner 
remains stationary (ignoring any velocity component perpendic- 



ular to the line-of-sight). The three-dimensional kinematic mo- 
tions can not be fully reconstructed with these data. 

5.6. The interstellar radiation field strength 

In order to estimate the effective interstellar radiation field for 
each interstellar cloud probed by the Upp Sco stars we con- 
structed a simplified model of a sheet of dust irradiated by sev- 
eral OB-type stars. The thickness of the dust sheet will roughly 
depend on the volume and column density of Hi. For Hh = 
100 cm" 3 and N(H) = 5 x 10 21 cirT 2 (E (B _ y) ~l mag) the thick- 
ness is ~16 pc. The dust sheet can be represented by a homoge- 
neous thin slab at a distance of 120 pc and a thickness of 20 pc. 
In this way, the distribution of individual clouds can be repre- 
sented by a single sheet, which is a valid assumption because 
(1) Tuv > 1 so photons are scattered frequently enough to loose 
most of their directional memory and (2) the distribution of in- 
dividual clouds has a surface area covering factor larger than 
unity. Property (1) assures that the radiation field strength I\jy is 
the roughly isotropic flux that impinges on the individual clouds 
making up the sheet. Aspect (2) assures that each line of sight 
through the representative sheet has approximately the same to- 
tal extinction, relevant for the attenuation of /uv- This ensures 
that the radiative transfer problem to be solved is that for a slab 
geometry. As eight OB stars contribute over 90% of the ISRF 
in this region t hese are included as the only source of the ion- 
ising radiation dSuiafha et al. 1 12005b . These stars illuminate the 
interstellar cloud from behind. In this particular model one star, 
HD 143275 (B0.3IV), dominates the effective ISRF, even while 
£ Oph (HD 149757) has the earliest spectral type. 

The radiative transfer model dSDaanslll995l ISnaansI fl9%l) 
takes into account both absorption and scattering. The effective 
optical depth Ty is computed from the observed E(b-v) and sub- 
sequently e~ T ' is multiplied by the individual stellar fluxes for an 
appropriate extinction curve for standard Milky Way dust (with 
Ry taken either as 3.1 or 4). The latter case is also considered 
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Fig. 15. W, 1 (5797)/W /i (5780) against W(CH)/E (B -v) (top) and 
W(CH + )fE(B-v) (bottom). Tentatively a positive trend can be dis- 
cerned between the DIB ratio and W(CH)/E(b-v) (top) but not 
between the DIB ratio and W(CH + )/E (B -v) (bottom). 
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Fig. 16. W /1 (5797)/W /1 (5780) vs. W(CH)/W(CH + ). There is a 
tentative trend for decreasing DIB ratio with increasing CH + /CH 
ratio. 



since the Upp Sco region contains sightlines with high Ry values 
for the dust extinction (see Sect. 15 .81 . In addition, this method is 
also applied to compute the amount of back-scattered radiation. 



For a sheet geometry, it is possible to express 7uv (in units of 
the Draine field) as a function of Ry and cloud position R. 



Ajv = 4.7[(/?+a)/a]"exp 



-6.9 



Z(B-V) 

0.20 



(fr)<» 



if) ■ 10 



-2 



(1) 



In this parametrisation of the radiative transfer grid the param- 
eter a depends on the distance to HD 143275 (a = 3.4 for a 
distance of 123 pc). E(b-v) is the individual extinction of the 
cloud. The cloud position R is set between and 20 pc (0 pc be- 
ing the sheet edge closest to the observer). Hence, /uv ~1 at the 
shielded edge, close to the mean Galactic value, and /tjv ~200 at 
the bright edge. This latter value is relatively high and depends 
on the distance to HD 143275. Placing this star 1 pc further away 
results in a = 5.1 and /tjv decreases by a factor two. With the 
nominal values, /tjv < 20 for about one third of the sheet struc- 
ture. Eq. Q] allows a range of impinging field strengths which 
has subsequently been used to set up a grid of chemical m od- 
els (including non-thermal production of CH; IS p aansl l 19951) for 
a given measured extinction, to determine which model clouds 
yield the best match to the available data. Thus effectively, for 
a given E^b-v) we extracted the /uv reproducing best the ob- 
served CH and CN, where the derived /tjv is also constrained 
by the observed upper limits for CN. In this, CH + has been ex- 
cluded because it is well known that canonical chemical models 
under produce its abundance by about two orders of magnitude. 
Turbulent dissipation and/or shocks are likely needed in the (en- 
dothermic) formation of CH + . To first order, /tjv/»h is the con- 
trolling parameter for the chemical and thermal balance. So an 
increase in density by a factor of 2 corresponds to an increase 
in /uv by a factor of 2. Due to the limited information avail- 
able for each line-of-sight we adopted a generic density = 
300 cirT 3 . This is representative of a cloud that is slightly denser 
than the ambient medium, the dust sheet, in which it is embed- 
ded. In other words, the clouds do not fill the region and the 
sheet is seen as a patchy complex of individual, but connected, 
clouds scattered in distance. Previous detailed modeling of the 
Upp Sco line-of-sight towards HD 147889 shows that this is 
likely a conservative lower limit for sightlines probing the denser 
parts of the p Oph cloud. With detailed model ing, including ob- 
servati onal constraints for additional species, iRuiterkamp et al.l 
(12005b found a density of 1200 cm 3 and an /uv ~10 for this 
line-of-sight. Increasing the input density by a factor of four in 
the model above for HD 147899 would give a revised /uv of 
6, already in better agreem ent with the detailed analysis. Also, 
IZsarpo & Federmanl (12003b found C i densities between 100 and 
300 cirT 3 for HD 143275 and HD 147165, but lower values, -50 
to -200, for HD 144470 and HD 144217. However, as C i traces 
the purely atomic phase of clouds, it is likely that this yields 
lower densities than for the molecular/shielded parts as traced 
by e.g. C2, CN and CH. Clearly, the simplifications introduced 
in the model presented in this work do not fully incorporate all 
the intricacies of a full-fledged analysis. However, the strength 
of this model, which relies only on the CH and CN abundance, 
is in giving statistically relevant predictions of the ISRF for a 
larger dataset for which only limited information is available. 
For accurate equivalent width measurements of both CH and CN 
the computed /uv has an uncertainty of approximately 25%, not 
including any unknown systematic effects. Uncertainties in the 
density, D = dn/n, propagate into /uv as D 1 ^ 2 . If only CH is de- 
tected the value of /uv should therefore be considered indicative 
only (like a model dependent lower limit). The resulting inter- 
stellar radiation field strengths are presented in Table [2] 
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For diffuse clouds the ISRF can also be estimated from 
steady-state gas ph ase chemistry (see e.g. IWeltv et all 120061 
iRitchev et a l. 2006): 



N(CH+) 
' N(CH) J - 



(2) 



which is valid for non-thermal CH production and for small val- 
ues of /h, . On the other hand, r otational excitation modeling o f 
H? gives (see e .g. lJural dl975h : iBlack & van Dishoeckl dl987h : 
iLee et al. (200^)): 



log(« H // uv ) oc log / 



(3) 



which is appropriate for n(H 2 ) << n(H) » rin (but the linear- 
ity holds also for higher N(H 2 ) {e.g. ILee et al.ll2007h . Note that 
Eqs. [2] and [3] show a n opposite dependence of /uv/"// on /h,. 

In addition, UV pumping can produce excited H* leading 
to an enhancement in the production of CH + via C + + H* 2 — » 
CH + + H. Therefore, we compare the independently obtained 
values for CH + and /uv, as well as CH/CH + and /uv to inves- 
tigate whether this process is important. Fig. [TBI illustrates that 
sightlines with high CH + abundances show only moderate val- 
ues for /uv (i.e. less than 10), while sightlines with high /uv 
larger than ~10) all show low-to-normal CH + abundances. From 
this relation it appears that a strong ISRF (/uv > 10) does not 
lead to enhanced CH + production, possibly because the molecu- 
lar hydrogen abundance of these sightlines is too low. However, 
W(CH + )/E(B-y) peaks at /uv = 5 which may reveal a delicate 
balance for the presence of UV pumping at intermediate /uv and 
moderate /h, . On the other hand, Fig.[19]shows that the CH/CH + 
ratio drops rapidly for /uv > 4. Thus despite a lower total CH + 
abundance (per unit reddening) for higher values of /uv the rela- 
tive production of CH + with respect to CH increases. This could 
be due to more efficient production of CH + or less efficient for- 
mation of CH in these low density, strongly UV exposed envi- 
ronments. The latter is indeed expected as N(CH) correlates with 
N(H 2 ) whose relative presence reduces also with increasing /uv 
(see below). UV pumping may thus contribute significantly to 
CH + formation only in diffuse clouds with sufficient abundance 
of both H 2 and UV photons. Other mechanisms, like turbulent 
dissipation of mechanical energy, could also be important for 
CH + formation in this region. 

Fig. [20] (top panel) reveals an evident inverse relation be- 
tween the strength of the ISRF, /uv, and the molecular hydro- 
gen fraction, /h, . This effect of lower /uv for interstellar clouds 
with higher molecular fractions (and thus more efficient shield- 
ing of the UV radiation) is expected from Eq. [3] This figure 
also illustrates the general trend that the cr-type clouds have 
a higher /uv and a lower molecular content /e,, while ^-type 
sightlines have a higher molecular content and are exposed to 
a weaker ISRF. Note that although a few sightlines with low 
/h, and higher /uv values were classified as £-type based on the 
observed W(5797)/W(5780) ratio, there is a clear separation - 
based on physical condition s - between the cr- and £-type sight- 
lines. IWeselak et all ( 120041) also show a similar distinction be- 
tween cr and £ at /h, ~ 0.4 (although their <x-£ classification 
is based on central depth ratios resulting in a slightly different 
division between the two types). 

The linear relation between log(/n,) and log(«////uv) 0-<? 
Eq. in Fig. [20l (bottom panel) can be compared directly to 
Fig. 2 in ILee et all (120071) showing indeed a close relation be- 
tween the molecular fraction and the ratio of hydrogen density 
over radiation field strength, nn/Iuv- This relation is sensitive 
to the total Ho column density but does not depend strongly on 




10 15 

ISRF [I_UV] 

Fig. 18. The CH + line strength per unit reddening is shown 
as a function of the ISRF strength, /uv- Highest values for 
/uv are found for low CH + abundances per unit reddening. 
Because N(CH) oc N(H 2 ) and E^b-v) 00 H 1 Eq. [2] gives /uv/"h 
oc N(CH + )/E(£_v')- There is some evidence for enhanced CH + 
production (UV pumping?) in clouds with moderate /uv ~ 5. 



the hydrogen particle density, nn- Nonetheless, knowledge of 
the latter value (either estimated or derived from complementary 
data) is required to derive /uv - Non-thermal H 2 excitation due to 
turbulence can mimic UV pumping and thu s alter the rela tion be- 
tween the model /uv and the observed /h, dSpaanslll995l) . Thus, 
in addition to deriving the effective ISRF strength in interstellar 
clouds (averaged along the line-of-sight) in Upp Sco with the 
model above we can use measurements of the /15780 DIB, CH 
and CN absorption line strengths to estimate N(H 2 ), N(Hi), and 
the molecular hydrogen fraction, /h, . 

Figure ED shows the dependence of the W(5797)/W(5780) 
ratio on the ISRF, /uv- In general, sightlines with low /uv val- 
ues are £-type sightlines for which also CN has been detected. 
And lines-of-sight for which we find high values of /uv have, on 
average, lower values for the W(5797)/W(5780) ratio. The few (- 
type clouds with high /uv have in fact DIB ratios that are close to 
the average ratio used to discriminate between cr and £-type en- 
vironments. On the other end, there are also a few cr-type clouds 
associated with a weak ISRF. We recall that a change in the den- 
sity will give an equal change of /uv, which could consequently 
shift individual sightlines to either lower or higher /uv, thus in- 
troducing additional scatter. Fig.|2T|is consistent with a 5797 A 
DIB carrier which requires sufficient protection from UV radia- 
tion in order to survive in the diffuse ISM, while conversely the 
5780 A DIB carrier requires UV photons for excitation (possi- 
bly because the carrier needs to be ionised in order to absorb at 
5780 A). At this point it is important to note that the sightlines 
with higher inferred /uv all rely on CH measurements only and 
should therefore be considered indicative. Also, the average /uv 
values for respectively cr and ( type sightlines are within lcr of 
each other (where the mean of /uvo- is two times the mean of 
/uvf)- Higher sensitivity data of CN transitions in Upp Sco are 
required to accurately probe /uv throughout the region. In that 
case, subsequent comparisons with accurate CH + line-widths (to 
determine the Doppler velocity parameter b) could be used to 
distinguish between the production of CH + in shocks {c.q. tur- 
bulent media) and the effect of UV pumping on enhanced abun- 
dances of CH + (see also Sect. 5.3). 
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Fig. 19. The CH over CH + line strength ratio is plotted as a func- 
tion of the ISRF strength, /uv- There is a drop in this ratio (i.e. 
enhanced CH + or reduced CH production) for stronger radiation 
fields (/uv)- This trend is consistent with non-thermal production 
of CH, otherwise no trend would be expected. 



5.7. Spatial distribution of DIBs and DIB ratios 

Our dataset provides a unique opportunity to investigate the 
scatter on the linear relation between DIB strength and redden- 
ing by dust, in particular with respect to its spatial distribution. 
Therefore, the equivale nt width per unit red dening is plotted on 
the infrared dust map dSchlegel et al.l[l998t Fig. l22l i. We show 
only the results for the 5780 and 5797 A DIBs. The equivalent 
widths for the 6196, 6379, and 6613 A DIBs behave similarly to 
the 5797 A DIB, but due to the larger relative uncertainties in the 
measured equivalent widths are not discussed further. 

From Fig. [22] it is apparent that the spatial behaviour for 
the two DIBs is different. In order to visualize this effect we 
show the spatial distribution of the W^(5797)/W,i(5780) ratio in 
Fig. [23] It can be seen that cr-type sightlines are more frequently 
probed towards the region westward of the p Oph cloud which 
has a low dust content, while predominantly ^-types are ob- 
served to wards the high dust c olumn density pOphiuchus cloud 
complex. ISuiatha et alj ( 120051) showed that 5 of the 8 stars pro- 
ducing 90% of the local ISRF are situated in this low-dust region 
at distances between 109 and 141 pc. Here, the strong stellar 
winds have blown out most of the dust and are now impinging 
on the west side of the p Oph cloud. 

To study this difference we selected two regions, a region 
free of dust emission and one with strong dust emission, respec- 
tively. The first region is centred on the p Oph cloud (showing 
a high dust column), while the second region is centred on the 
region west of p Oph scarce in dust emission (Fig. [23] black and 
white squares, respectively). In line with the observed dust den- 
sity, UV field strength and molecular H2 fraction, 80% of the 
lines-of-sight in the selected "high-density" (i.e higher dust col- 
umn and higher /h,) region are classified as £ and 95% of the 
sightlines in the "low-density" (i.e low dust column density and 
low fn 2 ) region are designated cr-type. 

The weighted mean and the associated error of the DIB ra- 
tio, Iuy, and /hi are calculated for the "low-density" and "high- 
density" selected regions, as well as for the total dataset (all 
sightlines with available values) and are given in Table [6] The 
two regions differ significantly from each other, and both regions 
show deviations from the overall weighted mean. For the "low- 
density" region the DIB ratios peak at about 0.20 and with a 
distribution width of about 0.05, while the "high-density" DIB 
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Fig. 20. (top) The model /uv is plotted as a function of inferred 
/h 2 ■ (bottom) Alternatively, the molecular hydrogen fraction can 
be plotted as a function of the ratio of the hydrogen density of 
the UV radiation field strength, for our general value of vih - 
300.0 cm -3 . As expected from iBlack & van Dishoeckl ( 1 19871) . 
log(/H,) is directly proportional to log(n#/7uvX where the in- 
tercept of this rel ation depends on the total H2 column density 
(see also Fig. 2 in lLee et ail d2007h ). Ranges for inferred N(Ht) 
are indicated by different symbols. The horizontal arrow gives 
the change in log(n#/7uv) for an increase or decrease of iih by 
a factor of 2. The vertical arrows on the data points indicate the 
correction of the inferred fraction to the directly observed frac- 
tion. Note that the inferred molecular fraction, fy 2 , directly de- 
pends (non-linearly) on the ratio of W(CH) over W(5780), and 
the /uv depends also on CH, as well as CN and Ecs_y). 



ratio distribution peaks at about 0.45 with a wider width of about 
0.10 (Fig. |23). The mean value of 0.26 ± 0.01 for the DIB ratio 
has been adopted to make the distinction between cr and type 
clouds. Fig. [20] shows that the mean molecular hydrogen frac- 
tion, /h 2 = 0.34 ± 0.21 (Tabled, provides an alternate - com- 
plementary - way to distinguish between cr and £ sightlines (e.g. 
Fig. [20] top panel). 
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Fig. 21. The W(5797)/W(5780) DIB ratio is plotted against the 
ISRF, Tuv- High DIB ratios (i.e. () correspond to a lower Zyv- 
cr-type clouds show a similar range in 7uv- 

Table 6. The weighted mean and associated error of 
W(5797)/W(5780), 7 UV , and f Hl for all sightlines (for which val- 
ues are available) and for the "low-density" and "high-density" 
regions. 



weighted mean 


all si, 


jhtlines 


"high-density" 


"low-density" 


W(5797)/tV(5780) 


0.26 


±0.01 


0.36 ± 0.03 


0.16 ± 0.01 


/uv 


6.4 ± 


2.5 


4.2 ± 2.5 


8.1 ± 6.2 


/h 2 


0.34 


±0.2 


0.44 ± 0.2 


0.15 ± 0.1 



5.8. R v andW A (5797)/W A (57B,0) 

For Milky Way lines-of-sight a weak relation exists between 
the total -to-selective visual extinction ratio, Ry, and the UV ex- 
tinction (ICardelii et al.ll989tlFitzpatrick & Massal2007h . For in- 
creasing Ry the far-UV absorption decreases, suggesting that 
fewer small dust particles / large molecules absorb in the 
far-UV, and implying a dust size distribution shifted towards 
larger grains. Therefore, Ry is sometimes used as a tracer 
of high density ISM, where grain growth is more significant; 
Ry - 3.1 for the average diffuse ISM, wh ile Ry = 4 - 6 
(ICardelii & Wallersteinl fl986l ICardellil fl988h for dense ISM. 
IWhittetl (11971 " showed that sightlines penetrating the pOph 
cloud reveal a higher than average Ry. Extinction curves with 
a higher Ry have less steep far-UV rise, which is associated to 
a lack of small particles which is often attributed to enhanced 
grain growth in denser clouds. However, an obvious correlation 
between R v and n H , N(CH), N(CN) could not be identified. The 
observed differences in Ry towards Upp Sco showed that the 
dust size distribution is not homogeneous throughout the primor- 
dial cloud that formed the association, or that the dust has been 
processed differently (for example, due to destruction in UV ex- 
posed regions, or grain growth in denser areas) in the various 
parts of the association. 

In an attempt to reveal any dependence between the UV 
radiation field and the dust grain size distribution we plot 
W / j(5797)/W / i(5780) versus Ry for 23 targets in our sample 
(see Table [TJ. Large error bars on Ry and the large differences 
between values obtained by different authors prevent us from 
discerning any significant trends between the DIB ratio and 




Fig. 22. Circles correspond to W / j(5780)/E(b_v) (top) and 
W A {5191)[E (B -V) (bottom). Background: IRAS 100 fim dust 
map. Targets are indicated by their target number as defined in 
Table Q] White circles represent cr-type sightlines and black cir- 
cles f-type sightlines, except those labeled "i" which are of in- 
termediate type. Targets labeled with "*" have relative errors on 
the ratio between 30 and 60%, while errors larger than 60% are 
omitted. The circle sizes for W(5797)/E(#_y) are multiplied by a 
factor 2 with respect to those for W(5780)/E(£_y). 



Ry. Th is is consistent with the earlier work of ISnow & Cohenl 
(fl97l . 



6. Conclusions and summary 

Lines of sight can be designated via W(5797)/W(5780) as ei- 
ther ^-type (sightlines penetrating cloud cores) or cr-type (sight- 
lines probing cloud edges) (or intermediate). We investigated 
the spatial variation of the DIB strengths W(DIB) and the 
W(5797)/W(5780) DIB ratio, and their dependence on redden- 
ing in 89 lines-of-sight within a field of 20°x20° probing the 
small scale variations in the gas and dust in the Upp Sco associ- 
ation. This represents an in-depth multi-object study of a well- 
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149438 



Fig. 23. The circle sizes is proportional to W, i (5797)/W,,(5780). 
Black and white circles indicate £ and cr sightlines, respectively. 
Seven of the eight stars (one is outside the map) generating 90% 
of the local 1SRF (ISujafha et al.ll2005h are indicated by bold cir- 
cles with italic numbers. The "high-density" (pOph cloud) and 
"low-density" regions selected for comparison are delineated by 
the black and white square, respectively. Intermediate classified 
lines-of-sight are labeled with an "i". Note that for HD 149438 
no DIBs were detected, but it is included here for its large con- 
tribution to the local ISRF. Targets labeled with "*" have relative 
errors on the ratio between 30 and 60%, while errors larger than 
60% are omitted. All other targets have errors smaller than 30% 
(see also Table©. 



□ all sightlines 
1 "low density" 
"high density" 



_hL± 



0.1 0.2 0.3 0.4 0.5 

W(5797)/W(5780) 

Fig. 24. Histogram of the W(5797)/W(5780) ratio for the lines- 
of-sight in the low-density (blue) and high-density (red) region, 
and for all sightlines in the sample. The vertical dashed line sets 
the division between cr and £-type cloud, (see electronic version 
for colour version). 



studied interstellar cloud complex, providing a valuable statisti- 
cal dataset. These data cover a wide range in dust column den- 
sities (from zero up to four magnitudes of visual extinction) to 
track the sensitivity of the DIB carrier molecules in relation to 
their local environment. 

Our results provide evidence that on average the DIB 
strengths in Upper Scorpius are linearly proportional to the 
reddening, closely following the general relation observed for 
the Galactic diffuse ISM. In addition, we showed that the 
scatter on these relationships, expressed for example via the 
W(5797)/W(5780) and W(DIB)/E (B _ V) ratios, is significant and 



Fig. 25. W(5797)/W(5780) is plotted against the total-to- 
selective visual extinction ratio Ry for 23 stars. The low and high 
density sightlines are in dicated by squa res and circles, respec- 
tively. Values are from Wegnei J 20021) [W0 21 ; IValencic et alJ 
(l2004h [V04], and lFitzpatrick & Massal d2007l) [FM07]. 



can be attributed to variations in the local physical conditions, in 
particular the interstellar density and the radiation field. 

We found that making a distinction between cr and f type 
sightlines clearly improved the relation between 5780 A DIB 
strength and the amount of dust, E(b-v)- The improvement for 
the other DIBs is less pronounced although still significant, sug- 
gesting that particularly the 5780 A DIB is sensitive to variations 
in the local conditions of the interstellar gas and dust. 

The CH and CH + molecules are detected in 53 out of 89 
sightlines, whereas CN and Cai are detected in 15 and 31 sight- 
lines, respectively. CH traces mainly cold (UV shielded) cloud 
material in the line-of-sight which is confirmed by the improved 
regression coefficient with respect to E^b-v) for £-type sightlines 
(r = 0.83) compared to its cr-type equivalent (r = 0.53). In line 
with previous results we observe no significant velocity differ- 
ence between the CH and CH + lines in the Upp Sco region, 
thus finding no evid ence for the production of CH + in shocks 
dRitchev et al.ll2011l) . 

CN traces the dense cloud cores and is only detected in 
type sightlines. The presence of (presumably very) low column 
densities of CN in cr-type sightlines remains to be confirmed. 

The VK(5797)/VK(5780) DIB ratio is more strongly correlated 
with W(CH)/E (B _v) than with W(CU + )/E (B -v)- This result con- 
firms that f-type clouds (i.e. dense and A5191 favoured) have a 
higher molecular content (e.g. higher abundances of CH and H2). 
The CH/CH + ratio is used to discriminate between quiescent and 
turbulent regions, confirming the division between cr and ^-type 
environments. This DIB ratio is thus related to the abundance of 
molecular species and dust properties inside the cloud cores. 

We found a significant difference between the mean 
VK(5797)/VK(5780) ratio towards the dense pOph cloud and 
towards the lower dust density region to the west. The high 
dust density region consists predominantly of f-type sightlines, 
whereas (as expected) the low-density region contains mainly cr- 
type lines of sight. This distinction was used to assign all sight- 
lines to either cr or £ type. 

Radial velocities of atomic and molecular absorption lines 
reveal one main interstellar component (with possibly unre- 
solved small-scale structure) at approximately -9 km s _1 , corre- 
sponding to the Upp Sco dust sheet at a distance of 1 10 to 150 pc. 
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Another, much weaker, velocity component at «-22 km s _1 , is 
linked to the tenuous foreground dust sheet at 50 pc. 

Asymmetries and widths of the Na i and K i line profiles sug- 
gests that the main velocity component may be resolved in nar- 
rower components at higher spectral resolution. The radial ve- 
locities for individual sightlines reveal velocity gradient in the 
dust sheet associated to Upp Sco. This could be seen as a differ- 
entially moving sheet, with the south-west corner fixed and the 
north-east corner moving towards the observer. For the sightlines 
with the strongest DIBs, the atomic, molecular, and DIB radial 
velocities are similar within uncertainties, indicating that these 
species spatially co-exist with each other. 

A simplified radiative transfer model of the main Upp Sco 
dust sheet was used to compute the effective interstellar radiation 
field (ISRF) from CN and CH observations. Five (of the eight) 
stars that strongly contribute to the local ISRF are located in the 
lower-density region of the Upp Sco cloud complex, substantiat- 
ing the inferred strong radiation field. The most accurate values 
for 7tjv are derived if both CN and CH are detected, and there- 
fore these sightlines correspond to those probing denser f-type 
regions. From the observed peak in the CH + /E^_y) ratio for 7uv 
ss 5 we can infer that UV pumping plays a role in the production 
of CH + in the Upp Sco region for diffuse clouds with sufficiently 
strong UV field and significant abundance of H2. The drop in 
CH/CH + with increasing 7uv is probably due to lower CH abun- 
dances for the more tenuous, UV exposed clouds. In fact the 
effective radiation field, 7uv, is an inverse function of the molec- 
ular hydrogen fraction, /h 2 , derived from N(CH) and W(5780). 
And, log(/Hi) is, as expected, proportional to log(/i#//ijv)- The 
W(5797)/ W(5780) ratio versus 7tjv relation in Upp Sco confirms 
that most of the sightlines with low /tjv values are of f-type 
for which also CN is detected, while only few f-type sightlines 
have large 7uv- The average Iyy for cr-type sightlines is higher 
than for ^"-types, but the standard deviation is relatively large for 
the poorly determined cr-type sightlines due to limited CN data. 
Exceptionally high S/N spectra are needed to probe the CN con- 
tent of these most tenuous diffuse clouds. 

Our results indicate that the relative abundance of the /15780 
DIB carrier increases with respect to the /I5797 DIB with an 
increase in the effective strength of the interstellar UV radia- 
tion field. A possible scenario is that a sufficiently large UV flux 
ionises the molecule whose cation has an electronic transition at 
5780 A. The /15797 DIB carrier, on the other hand, is more ef- 
ficiently destroyed in regions with high 7uv and is only able to 
survive in the deeper layers of the cloud. The behaviour of the 
/I/I5780 and 5797 DIBs is thus strongly dependent on the local 
effective UV radiation field (and consequently /h,) as is to be 
expected for molecular carriers. 

The scatter on the relation between DIB strength and red- 
dening can ultimately be translated to monitor the life cycle of 
the carriers as they get formed, excited, and destroyed through- 
out the complex interstellar cloud structure associated with the 
p Oph cloud and the Upp Sco association. The main driver for 
this cycle is the UV radiation field. More information is needed 
on the exact conditions of these diffuse clouds to disentangle ad- 
ditional (secondary) drivers. Measurements of additional atomic 
and molecular lines, such as C2, CN, H2, and CO, can improve 
the constraints on the chemical networks and physical conditions 
of individual clouds and thus improve our insight in the effects of 
these parameters on the abundance of DIB carriers. In addition 
to the two strong narrow DIBs at 5780 and 5797 A discussed 
in this work, accurate - high S/N - observations of additional 
narrow and broad diffuse bands are essential to fully probe the 
effects of changing environmental conditions on the ensemble 



of carriers giving rise to the DIB spectrum. Moreover, additional 
detailed studies of the Upp Sco region need to be complemented 
with similar studies of other specific Galactic and extra-galactic 
regions to probe a variety and range of different conditions that 
affect formation and destruction of the DIB carriers. 
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Appendix A: Error on equivalent width 

For the equivalent width measurement (W = ( Fc F F '' )), 
two integration borders A\ and A2 are selected for which the av- 
erage intensity in the continuum is derived in a window of width 
AA C = 1.4 A centred on the borders (Fig. IA.lt . 

The broader DIBs at 5778 and 5796 A underlying the strong 
5780 and 5797 A DIBs are largely removed from the spectra in 
the local continuum normalisation. They are very weak, even for 
high E(B_y), and do not significantly contaminate the W mea- 
surements of the 5780 and 5797 A DIBs. 



0.95 
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HD146029 
X D | B - 5780 A 
E B _ V = 0.15 mag 
S/N = 235.1 

= 96.7 ± 3.3 mA 



5774 5776 5778 5780 5782 5784 5786 

Fig.A.l. Typical example illustrating the determination of the 
DIB equivalent width. 



Taking into account both the noise of the measured flux 
(c"noise) an d the uncertainty in estimating the continuum (<x c ), 
the error on the equivalen t width is given according to 
IChalabaev & Maillardl(ll983l) [CM83], in our notation: 



2 _ 2 , 2 

^loise ^"continuum 

via equation A.9 from CM83 : 
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(A.l) 

(A.2) 
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SNR 

where we used the integration range A2 - A\ = M 5 A — AA\\ nf ,, 



SNR is the signal-to-noise ratio (per pixel), -0- 



AA/5A, 6 A is the spectral dispersion in A/pix, M the number of 
pixels in the integration range, and we take AA C « AA\i Be = AA. 
In Eq. I A. 61 the last term is approximately 2 as for shallow- weak 
lines W <sc AA = 3 A. We have 5A = 0.03 A, FWHM * 1 A (typ- 
ical for a narrow DIB), so that AA « 3 FWHM. This error esti- 
mate yields similar re sults as when applying the error estimate in 
the weak line limit bv lVollmann & Eversberd d2006l) (Eq. A.X): 

°> - V2 „.,„ - ™ VAAdA. 



For the strong narrow K 1 line the approximation in Eq. 
is not valid and the given cr w will be an upper limit to the true 
cr w . 

Appendix B: Contamination with the photospheric 
spectrum 

In order to evaluate possible contamination from stellar lines 
to the measured equivalent widths of the interstellar lines (both 
DIBs and molecules) we have checked synthetic spectra of B 
main-sequence [A0V, B5V, B2V, B0V, 09V] stars for the pres- 
ence ofatm^srAericJmes^^ of interest in this study 
(see iGummersbach & Kauferlll996l for details on the models). 
A brief overview of the expected contamination is given in 
TablelBTl 

For the 5797, 6 1 96 A DIBs no contamination above the 2 mA 
level is expected for these B-type main-sequence stars. The Nn 
line at 6610.6 A present in spectra of 09V to B2V stars is clearly 
separated from the 6613 A DIB at the spectral resolving power 
of FEROS. Both the On line at 4303.6 A and the Fen line at 
4303.2 A are also well resolved from the CH 4300.3 A line (the 
rest wavelength velocity difference is more than 200 km s _1 ). For 
CH (4300.3 A) the contamination in A0V spectra may be severe 
due to the presence of Ti n, however in that case also the Fe 11 line 
at 4303.2 A should be clearly present (with W ~2QQ mA). Strong 
Fe 1 lines appear in A0V spectra at 3872.5 A and 3878.6 A, how- 
ever, these are also clearly separated from the interstellar CN line 
(3874.6 A). The velocity profile plots cover a region of +0.6 A 
around the CN line. Also the weak to strong (W tot ai ~ 34 to 
~ 185 mA) C 11 lines around 3876.0 - 3876.7 A fall outside the 
plotted window. Only the Ne 11 line at 3875.3 A is relatively close 
to the CN line (with rest wavelength difference of 54 km s _1 ). 

There are several ways to investigate the potential prob- 
lem of contamination. For example, the line widths of all stel- 
lar lines are set by the rotational velocity (v sin i), whereas the 
di-atomic lines have intrinsically small widths, and can thus be 
distinguished from each other. Also, the radial velocities of the 
stellar and interstellar lines are (in virtually all cases) different. 
Therefore, any (di)atomic line with a velocity significantly dif- 
ferent from e.g. that of the Ki doublet (no contamination ex- 
pected) is suspect. Furthermore, in the case of relatively strong 
CH + line, the CH + transition at 3957.7 A can be used to ver- 
ify and confirm the abundance of CH + . For the relatively broad 
DIBs this differentiation due to velocity differences is less se- 
cure. Fortunately, the contamination expected at the DIB wave- 
lengths is very low. 

Table B.l. Overview of possible significant contamination from 
stellar lines to the interstellar line equivalent width measure- 
ments 

measured line contaminating line Spectral types contamination W 



5780 A 



Fen (5780.1 A) 



He 1 (5780.5 A) 
6379 A Nett (6379.6 A) 

CH (4300.3 A) Ti 11 (4300.0 A) 

CH+ (4232.5 A) Fe 11 (4232.8 A) 



Sin (4233.2 A) 



B5V-A0V 

<B5V 

09V 

B0V, B1V, B2V 

B5V 

A0V 

B2V 

B5V 

A0V 

B2V 



13 mA 

< 4 mA 
~2mA 

8, 22, 8 mA 
4 mA 
161 mA 

< 12 mA 
-98 mA 
-300 mA 

< 12 mA 
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Appendix C: Correlation plots 
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Fig. C.l. Equivalent width per unit reddening versus E^b-v) for 
the 5780 and 5797 A DIBs. Reduced chi-squared for linear 
fits are indicated in the respective panels. This plot reveals an 
intrinsic scatter on the relationship between the DIB strength 
and amount of dust in Upper Scorpius. Sightlines with E^b-v) 
< 0. 1 mag (orange/grey squares) are omitted for calculating the 
correlation coefficient; rsj^Q = -0.32 and ^5797 = 0.07. 
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Fig. C.2. Equivalent widths per unit reddening versus E(b-v) for 
the 6196, 6379, and 6613 A DIBs. Reduced chi-square for lin- 
ear fits are indicated in the respective panels. This plot reveals 
an intrinsic scatter on the relationship between the DIB strength 
and amount of dust in Upper Scorpius. Sightlines with E^b-v) 
< 0.1 mag (orange/grey squares) are omitted for calculating 
the correlation coefficient; r6i% = -0.41, ^379 = -0.17, and 

^6613 = -0.12. 
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Fig. C.3. Equivalent widths for CH (top) and CH + (bottom) are 
plotted versus W(6613), the strength of the 6613 A DIB. See 
Sect. 1541 



T 
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Fig.C.4. The W(6196)/W(6613) ratio plotted against E (B V) . 
Contrary to the (5797)/VK(5780) ratio discussed in Sect. l5?7l this 
ratio is less sensitive to reddening. 
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Appendix D: Line profiles 
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Fig. D.l. Observed profiles for the interstellar lines of CH + , CH, and Ki in sightlines probing the Upp Sco region. HD numbers 
correspond to those given in Table Q] In the bottom panel the weaker Ki doublet line at 7699 A is over-plotted (thick solid line) on 
top of the stronger (telluric contaminated) K i doublet line at 7665 A. The Ca i and CN lines are included only if potentially detected. 
The CH and CH + lines are shown for all lines-of-sight. The stronger CH + line at 4232 A is over-plotted (thick solid line) on top of 
the 3957 A line. 
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Fig. D.l. Continued. 
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Fig. D.l. Continued. 
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Fig. D.l. Continued. 
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Fig. D.l. Continued. 
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Fig. D.l. Continued. 
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Fig. D.l. Continued. 
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Fig. D.l. Continued. 
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Fig. D.l. Continued. 
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Fig. D.l. Continued. 
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Fig. D.l. Continued. 
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Fig. D.l. Continued. 
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Fig. D.l. Continued. 
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Fig. D.l. Continued. 
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Fig. D.l. Continued. 
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Fig. 6. The equivalent widths of the 6196, 6379, and 
6613 A DIBs in lines-of-sight previously measured and 
included here are compared with the measurements ob- 
tained i n this work. Literatu re valu es are taken from 
iHerbid (Il993l): ISeab & Snowl (Il995h: IGalazutdinov etai" 



2001: Ifhorburnet all d2003l): IGalazutdinov et al.l d200- 
ISollerman et al.l (120051) : iMegier et al.l d2001l l2005h . Linear 
regressions with the corresponding correlation coefficients r are 
shown in each panel. Slopes and intercepts of these regressions 
are 1 .06 and -0.8 mA, 1 . 1 3 and -3 . 1 mA, and 0.995 and - 1 .9 mA, 
for the 6196, 6379, and 6613 A DIBs, respectively. 
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Table 1. List of the observed targets. Cols. 1 to 8 give the Henry-Draper (HD) number, Hipparcos (Hip) number, spectral type, right 
ascension & declination, colo ur (B— V), intrinsic stellar colour (B — V)q , and reddening E^b-v)- Spectral types are fr om the Michigan 
Spectral Catalog of HD stars ( lHoukll982l ; lHouk & Smith-Moorell988l) . Intrinsic colours, (B- V)o, are taken from lFitzgeraldllT970i 
B and V photometric magnitudes are taken from Tycho-2 catalog and converted to the Johnson colour index B - V (using Bessell 
2000 and Mamajek et al. 2002); subsequently Eb-v = (B - V) - (B - V)p. Typ i cal err or o n E(b-v, is ~0.03 mag (see Sect. [3]) . Col. 9 



lists the total-to-selective visua l extinction ratio R v f rom lFitzpatrick & Massal d2007l) [fl.lWencic et alJ d2004t) fvl. lWegnerl (|2002|) 



rwl. lClavton & Cardellil(ll988h : ICardelli et al.l(ll989h fcl.lLewis et al.l (l2005h M. and lLarson & Whitted ((20051) Q21 Col. 10 gives the 
parallax distance d (pc) from Hipparcos dPerrvman et al 1 119971) . 



Nr. 


HD 




Hip 


Spectral 

tVT")P 


Right 
Ascension 


Declination 


(B - V) 


(B - V)n 


f-'B— V 


Rv 






HD 


1 38503 


76161 




1 5-33-1 8 91 

1J.JJ. 1 (J. y L 




-0.048 


-0.22 


0.17 






2 


HD 


139094 


76473 


B8IV 


15:37:06.89 


-26:29:32.6 


0.082 


-0.10 


0.18 




289.0 


3 


HD 


139160 


76503 


B9IV 


15:37:28.50 


-26:16:47.5 


-0.025 


-0.07 


0.05 




184.2 


4 


HD 


139486 


76633 


B9V 


1 5-39-00 06 

. ^ J .\J\) .\J\J 


-19:43:57.2 


0.033 


-0.07 


0.10 




193.4 


5 


HD 


139518 


76666 


B9.5V 


1 5-39-21 37 


-23 09-00 8 


0.021 


-0.04 


0.06 




93.7 




HD 


140543 


77131 


B 1 lah/Th 


1 5-44-56 66 

L *J . I I . ^)\J .\J\J 


-21 -48-53 9 


-0.013 


-0.19 


0.18 


3.16' 1 ' 


735.3 


7 


HD 


141180 


77449 


B9III 


15:48:42.90 


-27-34-56 3 


0.002 


-0.08 


0.08 




442.5 


8 


HD 


141444 


77569 


B9.5V 


15:50:11.47 


-28:42:15.6 


0.12 


-0.04 


0.16 






9 


HD 


141774 


77677 


B9V 


15:51:29.84 


-20:35:14.5 


0.063 


-0.07 


0.13 




202.8 


10 


HD 


142096 


77811 


B3V 


15:53:20.06 


-20:10:01.3 


-0.032 


-0.20 


0.17 




109.3 


1 1 


HD 


142114 


77840 


B2.5Vn 


15:53:36.72 


-25:19:37.7 


-0.105 


-0.22 


0.11 




133.0 


12 


HD 


142165 


77858 


B5V 


15:53:53.91 


-24:31:59.4 


-0.023 


-0.16 


0.14 




127.1 


13 


HD 


142184 


77859 


B2V 


15:53:55.86 


-23:58:41.1 


-0.052 


-0.24 


0.19 




120.5 


14 


HD 


142250 


77900 


B7V 


15:54:30.11 


-27:20:19.1 


-0.070 


-0.13 


0.06 




162.9 


15 


HD 


142301 


77909 


B8III/IV 


15:54:39.53 


-25:14:37.5 


-0.074 


-0.10 


0.03 


3.39" 


139.7 


16 


HD 


142378 


77939 


B2/B3V 


15:55:00.36 


-19:22:58.5 


-0.030 


-0.22 


0.19 




191.2 


17 


HD 


142669 


78104 


B2IV-V 


15:56:53.08 


-29:12:50.7 


-0.192 


-0.24 


0.05 




125.5 


18 


HD 


142883 


78168 


B3V 


15:57:40.46 


-20:58:59.1 


0.002 


-0.20 


0.20 




139.7 


19 


HD 


142884 


78183 


B8/B9III 


15:57:48.80 


-23:31:38.3 


-0.018 


-0.09 


0.07 




122.2 


20 


HD 


142983 


78207 


B8Ia/Iab 


15:58:11.37 


-14:16:45.7 


-0.089 


-0.02 


-0.07 




157.2 


21 


HD 


142990 


78246 


B5V 


15:58:34.87 


-24:49:53.4 


-0.100 


-0.16 


0.06 


3.22"' 


149.7 


22 


HD 


143018 


78265 


B1V 


15:58:51.11 


-26:06:50.8 


-0.154 


-0.26 


0.11 




140.8 


23 


HD 


143275 


78401 


B0.3IV 


16:00:20.01 


-22:37:18.2 


-0.083 


-0.28 


0.20 


3.60'' 3.09" 2.80+0.12' 


123.2 


24 


HD 


143567 


78530 


B9V 


16:01:55.46 


-21:58:49.4 


0.078 


-0.07 


0.15 




131.9 


25 


HD 


143600 


78549 


B9.5V 


16:02:13.56 


-22:41:15.2 


0.085 


-0.04 


0.12 




121.1 


26 


HD 


143956 


78702 


B9V 


16:04:00.24 


-19:46:02.9 


0.140 


-0.07 


0.21 




241.5 


27 


HD 


144175 


78809 


B9V 


16:05:19.15 


-23:40:08.8 


0.053 


-0.07 


0.12 




131.2 


28 


HD 


144217 


78820 


B0.5V 


16:05:26.23 


-19:48:19.6 


-0.070 


-0.26 


0.19 


4.00' 2.72" 3.90+0.27' 


162.6 


29 


HD 


144218 


78821 


B2V 


16:05:26.55 


-19:48:06.7 


-0.007 


-0.24 


0.23 






30 


HD 


144334 


78877 


B8V 


16:06:06.38 


-23:36:22.7 


-0.077 


-0.11 


0.03 




149.3 


31 


HD 


144470 


78933 


B1V 


16:06:48.43 


-20:40:09.1 


-0.057 


-0.26 


0.20 


3.75' 3.35" 3.37±0.29 v 
3.71+0.y 3.20+0.13' 


129.9 


32 


HD 


144569 


78956 


B9.5V 


16:07:04.67 


-16:56:35.8 


0.160 


-0.04 


0.20 




170.9 


33 


HD 


144586 


78968 


B9V 


16:07:14.93 


-17:56:09.7 


0.102 


-0.07 


0.17 




159.5 


34 


HD 


144661 


79031 


B8IV/V 


16:07:51.89 


-24:27:44.5 


-0.052 


-0.10 


0.05 




117.6 


35 


HD 


144708 


79005 


B9V 


16:07:36.42 


-12:44:43.5 


-0.006 


-0.07 


0.06 




130.5 


36 


HD 


144844 


79098 


B9V 


16:08:43.72 


-23:41:07.5 


0.012 


-0.07 


0.08 


4.04" 


130.7 


37 


HD 


144987 


79199 


B8V 


16:09:52.59 


-33:32:44.9 


-0.086 


-0.11 


0.02 




132.8 


38 


HD 


145353 


79343 


B9V 


16:11:33.52 


-27:09:03.1 


0.012 


-0.07 


0.17 




177.0 


39 


HD 


145482 


79404 


B2V 


16:12:18.20 


-27:55:34.9 


-0.154 


-0.24 


0.09 




143.5 


40 


HD 


145483 


79399 


B9Vvar 


16:12:16.04 


-28:25:02.3 


-0.044 


-0.07 


0.03 




91.4 


41 


HD 


145502 


79374 


B2IV 


16:11:59.73 


-19:27:38.6 


0.009 


-0.24 


0.25 


4.10' 3.48"' 3.57+0.15' 


133.9 


42 


HD 


145554 


79410 


B9V 


16:12:21.83 


-19:34:44.6 


0.125 


-0.07 


0.19 


3.65+0.43 / 


130.2 


43 


HD 


145556 


79437 


B4II/III 


16:12:43.56 


-28:19:18.1 


0.015 


-0.15 


0.16 
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Table 1. continued. 



Nr. 


HD 




Hip 


Spectral 
type 


Right 
Ascension 


Declination 


(B-V) 
(mag) 


(B - V) 
(mag) 


Efl-y 
(mag) 


Rv 




(pc) 


44 


HD 


145631 


79439 


B9V 


16:12:44.10 


-19:30:10.4 


0.133 


-0.07 


0.20 


3.82±0.53' 2 




125.3 


45 


HD 


145657 


79473 


B9.5V 


16:13:09.51 


-29:38:50.1 


0.143 


-0.04 


0.18 






370.4 


46 


HD 


145792 


79530 


B61V 


16:13:45.50 


-24:25:19.5 


-0.005 


-0.14 


0.13 


3.72±0.75' 




138.9 


47 


HD 


145964 


79599 


B9V 


16:14:28.88 


-21:06:27.5 


0.001 


-0.07 


0.07 






105.8 


48 


HD 


146001 


79622 


B8V 


16:14:53.43 


-25:28:37.1 


0.022 


-0.11 


0.13 


3.14" 




141.6 


49 


HD 


146029 


79621 


B9V 


16:14:53.33 


-22:22:49.1 


0.060 


-0.07 


0.13 






238.1 


50 


HD 


146284 


79740 


B9III/IV 


16:16:26.69 


-24:16:55.0 


0.125 


-0.08 


0.20 


3.l¥ 




263.9 


51 


HD 


146285 


79739 


B8V 


16:16:25.17 


-24:59:19.5 


0.201 


-0.11 


0.31 


3.83±0.18' 




179.9 


52 


HD 


146331 


79771 


B9V 


16:16:50.63 


-25:51:46.7 


0.301 


-0.07 


0.37 






145.8 


53 


HD 


146332 


79775 


B3III 


16:16:52.55 


-29:44:37.4 


0.132 


-0.20 


0.33 








54 


HD 


146416 


79785 


B9V 


16:16:58.77 


-21:18:14.9 


0.006 


-0.07 


0.08 






125.2 


55 


HD 


146706 


79897 


B9V 


16:18:28.26 


-23:16:27.5 


0.103 


-0.07 


0.17 


3.36" 




170.6 


56 


HD 


147009 


80019 


RN 


16:20:03.98 


-20:02:41.5 


0.256 


-0.01 


0.27 






160.8 


57 


HD 


147010 


80024 


B911/III 


16:20:05.49 


-20:03:23.0 


0.126 


-0.08 


0.21 


2.81±0.45' 2 




143.3 


58 


HD 


147103 


80063 


B9/A0V 


16:20:30.26 


-20:07:03.9 


0.338 


-0.04 


0.38 








59 


HD 


147165 


80112 


B1HI 


16:21:11.32 


-25:35:34.1 


0.097 


-0.26 


0.36 


3.80' 3.22" 
3.60±0.14' 


3.86+0.52" 
4.22±0.20' 


225.2 


60 


HD 


147196 


80126 


B5/8V 


16:21:19.19 


-23:42:28.7 


0.127 


-0.14 


0.27 


4.03" 3.10' 




151.3 


61 


HD 


147648 


80338 


B8II 


16:24:02.89 


-25:24:54.0 


0.620 


-0.10 


0.72 


3.49" 




153.8 


62 


HD 


147683 


80405 


B4V 


16:24:43.72 


-34:53:37.5 


0.119 


-0.18 


0.30 






591.7 


63 


HD 


147701 


80371 


B5III 


16:24:21.32 


-25:01:31.4 


0.500 


-0.16 


0.66 


4.05' 3.60" 
4.03±0.09' 


3.86±0.24 v 


149.7 


64 


HD 


147888 


80461 


B3/B4V 


16:25:24.28 


-23:27:36.8 


0.253 


-0.19 


0.44 


4.03' 

4.08±0.13' 


3.89±0.20' 


136.4 


65 


HD 


147889 


80462 


B2HI/IV 


16:25:24.32 


-24:27:56.6 


0.750 


-0.24 


0.99 


4.20 c 3.68" 
4.19±0.073 


3.95+0.13' 

f 


135.9 


66 


HD 


147932 


80474 


B5V 


16:25:35.08 


-23:24:18.8 


0.254 


-0.16 


0.41 


4.58 4.58' 




128.9 


67 


HD 


147933 


80473 


BIV 


16:25:35.12 


-23:26:49.8 


0.170 


-0.2 


0.4 








68 


HD 


147955 


80493 


B9.5V 


16:25:50.66 


-26:34:06.3 


0.206 


-0.07 


0.28 






154.8 


69 


HD 


148041 


80552 


B91V 


16:26:46.59 


-33:04:23.1 


0.316 


-0.07 


0.39 






442.5 


70 


HD 


148184 


80569 


B2Vne 


16:27:01.43 


-18:27:22.5 


0.203 


-0.24 


0.44 


4.15' 




149.9 


71 


HD 


148499 


80778 


B9III 


16:29:36.73 


-27:34:13.7 


0.414 


-0.08 


0.49 






346.0 


72 


HD 


148579 


80804 


B9V 


16:29:59.20 


-25:08:52.0 


0.245 


-0.07 


0.31 


4.01±0.19' 




160.3 


73 


HD 


148594 


80819 


B8V 


16:30:15.68 


-27:54:58.5 


0.074 


-0.11 


0.18 


3.W 




166.9 


74 


HD 


148605 


80815 


B3V 


16:30:12.48 


-25:06:54.8 


-0.132 


-0.20 


0.07 


3.18+0.52 w 




120.5 


75 


HD 


148860 


80940 


B9.5III 


16:31:38.71 


-17:42:48.3 


0.132 


-0.05 


0.18 






207.0 


76 


HD 


149367 


81211 


B8/9IV/V 


16:35:15.21 


-26:28:42.1 


0.105 


-0.07 


0.18 


3.49±0.10" 






77 


HD 


149387 


81236 


B7II/III 


16:35:29.22 


-30:16:03.5 


0.141 


-0.12 


0.26 








78 


HD 


149438 


81266 


BOV 


16:35:52.95 


-28:12:57.7 


-0.198 


-0.30 


0.10 






131.8 


79 


HD 


149757 


81377 


09.5Vn 


16:37:09.54 


-10:34:01.5 


0.020 


-0.30 


0.32 


3.09' 

3.08±0.19' 


2.55±0.24' 


140.4 


80 


HD 


149883 


81487 


B9V 


16:38:36.29 


-26:59:21.6 


0.128 


-0.07 


0.20 






237.5 


81 


HD 


149914 


81474 


B9.5IV 


16:38:28.65 


-18:13:13.7 


0.231 


-0.04 


0.27 






165.0 


82 


HD 


150514 


81785 


B8HI 


16:42:17.41 


-20:11:13.3 


0.103 


-0.10 


0.20 








83 


HD 


150814 


81941 


B9.5V 


16:44:17.66 


-22:31:24.9 


0.129 


-0.04 


0.17 






183.2 


84 


HD 


151012 


82069 


B9.5V 


16:45:48.46 


-26:38:57.9 


0.044 


-0.04 


0.08 






110.1 


85 


HD 


151346 


82217 


B8II 


16:47:46.53 


-23:58:27.5 


0.386 


-0.10 


0.49 


3.70" 3.81±0.1tt f 


253.8 


86 


HD 


151496 


82271 


B9V 


16:48:36.25 


-21:56:28.1 


0.191 


-0.07 


0.26 






298.5 


87 


HD 


152516 


82734 


B2III 


16:54:38.40 


-21:52:46.6 


0.061 


-0.24 


0.30 








88 


HD 


152655 


82822 


B9111 


16:55:32.54 


-21:34:10.2 


0.100 


-0.08 


0.18 






425.5 


89 


HD 


152657 


82839 


B811 


16:55:44.38 


-25:32:02.1 


0.004 


-0.10 


0.10 






361.0 



HD 










Equivalent Width W (inA) 












Line strength ratios 


fs, 


Ajv 




subgroup 




.15780 


.15797 


,16196 


,16379 


,16613 


CH(4300) CH+(4232) 


CN(3875) 


Cai(4227) 


Ki(7699) 


Ki(7665) 


1 t"7CIT/ ) ^7QI1 


AO 1 yo//ioo 1 J 


CH/CH + 




138503 


98.0 ± 4.2 


24.4 ± 2.2 


15.2 ± 2.4 


15.9 ± 2.1 


35.9 ± 2.7 


0.2 


10.6 


_ 




48.4 ± 3.8 


63.6 ± 3.6 


0.25 ± 0.02 


0.42 ± 0.07 


0.02 


0.02 


(50) 


ml 




139094 


46.0 ± 3.0 


19.7 ± 2.2 


13.0 ± 1.9 


9.5 ± 1.9 


17.3 ± 1.7 


5.7 


4. 1 


5.0 


_ 


121.7 ± 3.9 


162.1 ± 3.4 


a a i i a az 

0.43 ± 0.06 


0.75 ±0.13 


1 .4 


0.53 


3.6 


£ 




139160 


_ 


10.8 ± 2.9 


<6.9 


9.3 ± 3.1 


21.2 ± 2.9 


1 .5 


7. 1 


_ 




25.6 ± 5.8 


39.5 ± 5.9 




< 0.33 


0.21 




3.5 






139486 


38.5 ± 6.9 


< 11.8 


< 11.8 


< 19.6 


< 15.3 






_ 




25.9 ± 9.5 


18.5 ± 8.3 


< 0.32 














139518 


< 17.6 


< 6 


<6 


< 5.1 


< 6 






_ 




<3.3 


_ 












tr 




140543 


123.0 ±5.1 


14.7 ± 4.0 


23.7 ± 3.7 


25.0 ± 3.4 


51.6 ± 4.4 






_ 




43.6 ± 6.8 


50.0 ± 8.0 


a i t i a m 

0.12 ± 0.03 


f\ A£. i A AO 

0.46 ± 0.08 








(X 




141180 


39.3 ± 4.6 


< 10.3 


11.7 ± 2.2 


< 12.6 


9.5 ± 2.6 






_ 




6.3 ± 5.1 


_ 


< 0.26 


Y.li ± 0.4U 












141444 


43.8 ± 6.3 


12.7 ± 3.7 


9.5 ± 3.5 


< 22.2 


14.7 ± 3.3 










49.4 ± 9.0 


51.8 ±4.8 


a ta . a aa 
0.29 ± 0.09 


f\ £. a i a ^o 
0.64 ± 0.28 








£ 




141774 


82.7 ±4.1 


15.6 ± 1.9 


10.2 ± 1.6 


9.7 ± 1.8 


TOO . T T 

28.8 ± 2.3 


1 .2 


2.5 


_ 




24.4 ±3.1 


22.7 ± 3.8 


A 1 A . A AT 
0.19 ± 0.02 


a o c i a r\£. 
0.35 ± 0.06 


0.48 


0.1 1 


1 0.j 


(T 




142096 


89.5 ± 3.1 


13.7 ± 1.6 


5.6 ± 1.2 


< 8.4 


22.9 ± 1.7 




0.45 


_ 


5.8 


23.0 ± 3.1 


19.0 ± 2.9 


A 1 C i A AO 

0.15 ± 0.02 


A . A AC 

0.25 ± 0.05 






(9.0) 


tr 


Low p 


142114 


65.5 ± 2.8 


12.5 ±2.4 


8.1 ± 1.2 


< 5.5 


20.9 ±1.4 


1 .3 






_ 


22.4 ± 2.6 


20.1 ±4.0 


A 1 A . A A A 

0.19 ± 0.04 


a oa i a 
0.39 ± 0.06 




0.15 


7.5 


<T 


Low p 


142165 


89.5 ± 3.4 


11.3 ± 1.6 


4.0 ± 1.2 


< 9.7 


14.3 ±1.4 






_ 




18.8 ±4.0 


15.7 ± 3.4 


AIT .A AO 

0.13 ± 0.02 


A 1 O i A AA 

o.zo ± o.oy 








{T 


Low p 


142184 


74.5 ± 2.7 


15.5 ± 1.5 


7.3 ± 1.2 


10.9 ± 1.8 


19.9 ±1.8 




2.0 


_ 




32.2 ± 3.2 


48.2 ± 4.4 


ATI , A AO 
0.21 ± U.U2 


A 1C i A A"7 

0. j6 ± 0.0 / 


1 .4 


A 1C 

0.25 


6.4 


iT 


Low p 


142250 


34.5 ± 4.5 


<8.5 


< 10.3 


< 11 


< 13.1 






_ 




8.7 ±4.2 


_ 


< 0.26 














142301 


61.0 ± 3.7 


14.1 ±2.8 


22.1 ± 2.4 


10.1 ± 2.2 


10.6 ± 2.6 


■ 








10.5 ± 2.4 


_ 


A T T . A AC 

0.23 ± 0.05 


^ AA i i \ ^ ! 

2.09 ± 0.56 








{T 


Low p 


142378 


87.4 ± 2.8 


14.7 ± 2.2 


10.7 ± 1.1 


6.9 ± 1.8 


17.4 ± 1.6 


0.9 


4.2 


_ 


_ 


22.6 ±2.1 


23.0 ± 3.7 


AIT . A AT 

0. 17 ± 0.03 


A tZ 1 i A AA 

o.6i ± o.oy 


0.21 


A AO 

0.08 


20.0 


(J 


Low p 


142669 


< 11.4 


< 5.5 


< 2.0 


<4.2 


< 3. 1 






_ 




4.5 ± 2.3 


_ 
















142883 


95.1 ± 2.2 


11.5 ± 1.4 


8.1 ± 1.1 


< 7.8 


24.4 ±1.7 


4.9 


10.9 


_ 




65.9 ±2.1 


84.1 ± 3.8 


AIT . A AT 

0.12 ± 0.02 


A 1 i A AC 

0.33 ± 0.05 


0.45 


0.32 


3.7 


<J 


Low p 


142884 


107.7 ± 3.1 


13.9 ± 1.8 


11.5 ± 1.3 


15.3 ± 1.2 


28.6 ± 2.2 


1 .7 


5.8 


_ 




35.3 ± 2.3 


44.8 ± 5.7 


AIT .A AT 

0.13 ± 0.02 


A A A i A A£T 

0.40 ± 0.06 


0.29 


0.13 


2.8 


tr 


Low p 


142983 


81.8 ±2.3 


9.7 ± 1.5 


8.4 ± 1.1 




27.0 ± 1 .5 






_ 




7.2 ± 2.6 


6.8 ± 3.3 


A 1 T I A AT 
0.12 ± U.02 


All i A f\A 

O.j 1 ± 0.04 








tr 




142990 


51.5 ± 2.0 


< 10.1 


12.9 ± 1.6 


< 6.9 


10.1 ± 1.5 


0.9 


4.6 


_ 




11.3 ± 2.0 


6.8 ± 2.9 


< 0.20 


1 "IO i A 1C 

l.zo ± 0.2j 


0.20 


A 

0.13 


j. 1 


tr 


Low p 


143018 


28.8 ± 3.0 


< 8.2 


<7.8 


< 11.5 


10.5 ± 2.0 






_ 


_ 


4.7 ± 2.9 


2.6 ± 2.4 


< 0.28 


A 1 A i A 1 1 

0. /4 ± 0.2 1 




[0.08] 








143275 


68.2 ± 3.1 


20.9 ± 2.2 


7.4 ± 1.5 


12.3 ± 1.5 


20.3 ±1.2 


1.3 


2.1 


_ 


_ 


30.9 ± 2.2 


41.3 ± 1.8 


A 1 1 i (\ f\A 

0.3 1 ± U.04 


A 1C i A AO 

0.J6 ± O.Oo 


0.62 


0. 14 [0.1 1 ] 


1 1.4 


r 

b 


Lowp 


143567 


162.5 ± 3.6 


19.8 ± 2.5 


12.1 ± 2.3 


13.2 ±2.0 


40.1 ± 2.7 






_ 




21.0 ± 3.1 


19.0 ± 3.6 


0.12 ± 0.02 


0.30 ± 0.06 








tr 


Low p 


143600 


81.3 ± 5.3 


17.7 ± 2.4 


8.6 ± 1.6 


13.8 ± 2.5 


19.2 ± 3.0 






- 


- 


30.0 ± 4.0 


46.2 ± 2.8 


0.22 ± 0.03 


0.45 ±0.11 








tr 


Low p 


143956 


178.0 ±4.5 


19.4 ± 2.8 


17.2 ±2.0 


17.8 ± 2.4 


49.9 ± 2.2 










23.1 ± 5.8 


29.5 ± 5.0 


0.11 ±0.02 


0.34 ± 0.04 








tr 


Low p 


144175 


35.8 ± 3.6 


<2.7 


11.9 ± 3.2 


11.4 ± 2.1 


< 5.7 










5.4 ±2.1 




< 0.08 


> 2.1 










Low p 


144217 


152.3 ± 2.1 


21.6 ± 1.7 


12.7 ±0.9 


15.6 ± 1.2 


44.4 ± 1.4 


1.4 


4.8 




1.2 


23.1 ± 2.0 


32.3 ±2.2 


0.14 ±0.01 


0.29 ± 0.02 


0.29 


0.08 


13.3 


o~ 


Low p 


144218 


158.8 ± 2.3 


17.6 ± 1.8 


14.0 ± 1.0 


19.2 ± 1.4 


42.9 ± 1.1 


0.8 


4.7 




1.5 


28.0 ± 2.0 


33.8 ± 1.5 


0.11 ±0.01 


0.33 ± 0.02 


0.17 


0.04 


20.0 


tr 


Low p 


144334 


52.9 ± 1.5 


12.0 ± 1.5 


18.5 ± 1.2 


8.2 ± 1.6 


10.3 ± 1.7 










4.5 ± 1.7 




0.23 ± 0.03 


1.79 ± 0.32 








CT 


Low p 


144470 


165.8 ±2.1 


28.0 ± 1.6 


19.0 ± 1.0 


26.2 ± 1.2 


58.2 ± 1.5 


3.0 


5.9 




1.5 


29.1 ± 1.4 


34.8 ± 2.0 


0.17 ±0.01 


0.33 ± 0.02 


0.51 


0.14 [0.13] 


7.0 


tr 


Low p 


144569 


86.7 ± 4.8 


11.5 ±2.7 


11.3 ±2.2 


18.4 ± 2.7 


25.5 ± 2.9 










29.7 ± 5.2 


30.0 ± 4.5 


0.13 ±0.03 


0.44 ±0.10 








tr 




144586 


81.7 ± 3.4 


15.0 ± 2.2 


10.0 ± 1.8 


11.1 ± 2.0 


30.8 ± 2.5 










27.2 ± 2.7 


39.7 ± 2.0 


0.18 ±0.03 


0.32 ± 0.06 








tr 




144661 


57.4 ± 1.6 


9.0 ± 1.4 


<5.6 


< 5.8 


6.9 ± 1.6 










4.1 ± 2.4 


2.4 ± 2.2 


0.16 ±0.02 


< 0.81 








tr 


Low p 


144708 


61.0 ± 2.6 


11.8 ±2.1 


4.8 ± 1.2 


8.5 ± 1.8 


17.3 ± 2.0 










13.1 ± 2.3 


17.5 ± 3.1 


0.19 ±0.04 


0.28 ± 0.08 








tr 




144844 




< 8.0 


<5.7 


< 8.6 


21.3 ± 1.3 










9.0 ± 1.8 


9.5 ± 3.0 




<0.27 










Low p 


144987 


< 7.0 


<7.8 


<4.6 


< 3.4 


7.6 ± 1.4 










4.7 ± 1.5 






0.61 ±0.13 












145353 


95.5 ± 2.9 


32.7 ± 1.7 


9.0 ± 1.4 


11.5 ± 1.7 


31.8 ± 1.9 


4.1 


7.1 




1.8 


48.9 ± 2.9 


57.0 ± 3.0 


0.34 ± 0.02 


0.28 ± 0.05 


0.58 


0.28 


5.1 


( 




145482 


< 10.3 


< 3 


< 1.4 


<4.4 


4.6 ± 1.6 










3.0 ± 1.8 






<0.3 












145483 


< 13.2 


< 12.2 


< 2.3 


<4.7 


4.9 ± 1.4 










7.5 ± 2.1 






<0.47 












145502 


166.2 ± 2.0 


35.0 ± 1.5 


15.8 ± 1.1 


33.6 ± 1.2 


57.0 ± 1.7 


2.0 


4.1 




3.8 


38.9 ±1.1 


55.5 ± 2.0 


0.21 ±0.01 


0.28 ± 0.02 


0.49 


0.10 [0.09] 


8.2 


tr 


Low p 


145554 


143.0 ± 3.8 


25.2 ± 2.0 


16.0 ± 2.0 


24.5 ± 1.8 


48.3 ± 2.0 


1.1 


2.7 




1.7 


39.2 ± 3.5 


47.9 ±4.1 


0.18 ±0.01 


0.33 ± 0.04 


0.41 


0.07 


17.4 


tr 


Low p 


145556 


108.5 ± 5.5 


24.0 ± 3.1 


14.9 ± 2.3 


11.3 ± 2.1 


36.9 ± 3.6 


1.8 


2.6 






27.9 ± 3.8 


17.6 ±4.0 


0.22 ± 0.03 


0.40 ± 0.07 


0.69 


0.13 


6.5 


cr 
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Equivalent Widlh W (mA) 












Line strength ratios 


/h 2 


hv 


a* 


subgroup 




,15780 


..15797 


,16196 


,16379 


,16613 


CH(4300) CH+(4232) 


CN(3875) 


Cai(4227) 


Ki(7699) 


Ki(7665) 


,15797/,15780 


,16196/,16613 


CH/CH+ 




145631 


129.9 ±3.7 


27.1 ± 2.1 


14.9 ± 1.7 


19.0 ± 2.0 


36.5 ± 2.3 


3.8 


3.8 


- 


2.9 


47.9 ± 2.3 


61.0 ±4.0 


0.21 ± 0.02 


0.41 ± 0.05 


1.0 


0.21 


5.4 






145657 


97.2 ±5.8 


34.7 ± 4.3 


10.6 ± 2.8 


29.4 ± 2.7 


57.0 ± 3.3 


5.4 


14. 1 


- 


5.7 


78.0 ±4.1 


93.7 ± 6.9 


0.36 ± 0.05 


0.19 ± 0.05 


0.38 


0.33 


3.4 


b 




145792 


66.0 ± 2.6 


12.8 ± 1.9 


8.9 ± 1.7 


14.1 ± 1.5 


14.6 ± 1.4 


1.5 


0.2 


- 


- 


14.4 ± 1.8 


12.9 ± 2.3 


0.19 ± 0.03 


0.61 ± 0.13 


7.5 


0.17 


7.7 




Low 


145964 


43.1 ± 2.3 


8.1 ± 1.4 


3.3 ± 1.1 


7.6 ± 1.3 


16.2 ± 1.9 






- 


- 


4.6 ± 1.7 


- 


0.19 ± 0.03 


0.21 ± 0.07 










Low 


146001 


68.8 ± 5.4 


10.1 ± 2.5 


11.7 ± 3.4 


< 10.3 


17.7 ± 3.3 






- 


- 


7.8 ± 3.5 


- 


0.15 ± 0.04 


0.66 ± 0.23 












146029 


96.8 ± 3.3 


19.9 ± 2.4 


6.5 ± 1.5 


< 11.3 


22.5 ± 1.8 






- 


- 


13.1 ± 3.0 


10.5 ± 3.5 


0.21 ± 0.03 


0.29 ± 0.07 










I 'iw'' 


146284 


99.7 ± 2.8 


22.8 ± 1.7 


9.6 ± 1.0 


<7.8 


28.8 ± 1.5 


5.S 


3.9 


- 


1.4 


56.3 ± 2.7 


58.2 ±2.6 


0.23 ± 0.02 


0.33 ± 0.04 


1.5 


0.34 






I iw^ 


146285 


158.5 ± 3.9 


23.1 ± 2.3 


9.9 ± 1.4 


13.2 ± 1.6 


43.1 ± 2.0 


5.7 


9. 1 


- 


- 


31.9 ±2.7 


34.1 ± 2.9 


0.15 ± 0.01 


0.23 ± 0.04 


0.63 


0.25 


3.3 




owp 


146331 


189.8 ± 3.4 


46.1 ± 2.5 


17.8 ± 1.5 


19.6 ± 1.3 


70.2 ± 3.3 


6.4 


3.0 


- 


3.9 


57.0 ± 2.7 


87.0 ± 3.0 


0.24 ± 0.01 


0.25 ± 0.02 


2.1 


0.24 


6.6 






146332 


103.0 ± 2.9 


65.3 ± 2.4 


15.3 ±2.2 


37.1 ±2.0 


69.2 ± 2. 1 


1 5 1 


3.8 


9.3 


3.0 


106.5 ± 4.2 


107.5 ± 2.4 


0.63 ± 0.03 


0.22 ± 0.03 


4.0 


0.57 


5.5 


r 

b 




146416 


63.9 ± 2.7 


10.4 ± 1.2 


4.9 ± 1.2 


< 3.4 


19.7 ± 1.8 






- 


- 


6.5 ±2.1 


3.8 ± 2.5 


0. 16 ± 0.02 


0.25 ± 0.06 






2.5 




owp 


146706 


45.2 ± 4.0 


27.7 ± 2.6 


4.0 ± 2.3 


<4.2 


17.2 ± 2.7 


7.2 


5.0 


- 


1.4 


53.9 ± 5.2 


49.1 ± 5.0 


0.61 ± 0.08 


0.23 ± 0.14 


1,4 


0.58 


2.4 


b 




147009 


123.7 ±4.3 


31.7 ± 3.5 


12.9 ±2.5 


23.5 ± 2.3 


51.7 ± 3.5 


8.0 


1 8.3 


- 


- 


47.6 ± 3.3 


52.4 ± 6.9 


0.26 ± 0.03 


0.25 ± 0.05 


0.44 


0.37 


4. 1 






147010 


- 


- 


- 


12.1 ±2.0 




7.0 


14.1 


1.9 


- 


50.5 ± 3.1 


55.8 ±5.2 






0.50 




3.4 






147103 


131.0 ±3.4 


71.3 ± 2.8 


16.1 ± 1.5 


28.8 ± 1.6 


56.7 ± 2.3 


23.8 


7.4 


8.6 


- 


75.0 ± 1.8 


78.4 ± 3.7 


0.54 ± 0.03 


0.28 ± 0.03 


3.2 


0.63 


] .6 


b 




147165 


234.5 ± 3.9 


55.2 ± 2.2 


18.6 ± 1.8 


25.7 ±2.1 




3.8 


4.8 


- 


- 


27.1 ± 2.9 


28.1 ±2.7 


0.24 ± 0.01 




0.79 


0.13 [0.05] 


9. 1 




p Oph 


147196 


31.0 ±3.0 


31.3 ± 2.4 


4.9 ± 1.7 


<9.9 


25.6 ± 2.0 


15.0 


4.9 


2.4 


2.1 


61.6 ± 3.7 


58.9 ± 3.6 


1.01 ± 0.12 


0.19 ± 0.07 


3.1 


0.80 


2.4 


r 

b 


p Oph 


147648 


253.9 ± 4.8 


113.8 ±3.4 


33.2 ±2.2 


72.0 ± 2.2 


142.6 ± 3.2 


2 1 .4 


7.4 


- 


- 


89.3 ± 2.3 


90.3 ± 6.1 


0.45 ± 0.02 


0.23 ± 0.02 


2.9 


0.45 


3.2 


b 


p Oph 


147683 


96.7 ± 3.2 


33.9 ± 1.8 


11.2 ± 1.7 


15.0 ± 1.6 


45.2 ± 2.1 


20.0 


1 8.0 


17.4 


2.4 


100.4 ± 2.2 


84.3 ± 4.7 


0.35 ± 0.02 


0.25 ± 0.04 


1.1 


0.65 


] .6 


b 




147701 


256.4 ± 4.7 


95.6 ± 2.3 


21.3 ±2.0 


51.3 ± 1.7 


65.6 ± 2.2 


27.4 


17 7 


19.3 


4.1 


102.1 ± 3.1 


97.1 ± 2.9 


0.37 ± 0.01 


0.33 ± 0.03 


1.6 


0.51 


2.4 


r 


p Oph 


147888 


196.9 ± 6.6 


49.7 ± 4.4 


19.2 ±4.0 


32.3 ± 2.8 


62.0 ± 4.3 


14.8 


6.9 


3.6 


6.1 


90.6 ± 7.0 


110.1 ± 9.9 


0.25 ± 0.02 


0.31 ± 0.07 


2.1 


0.41 [0.18] 


] .9 




p Oph 


147889 


343.7 ±4.1 


144.7 ± 2.9 


39.9 ± 2.8 


83.5 ±2.1 


165.4 ± 2.4 


45.5 


24.7 


25.6 


7.8 


83.0 ±2.1 


80.6 ± 3.0 


0.42 ± 0.01 


0.24 ± 0.02 


1.8 


0.56 


] .5 


b 


p Oph 


147932 


207.3 ±4.1 


60.3 ± 2.0 


14.4 ± 1.5 


25.8 ± 1.7 


57.1 ± 1.7 


1 3. 1 


4.9 


13.0 


3.6 


87.0 ± 3.4 


96.5 ± 2.5 


0.29 ± 0.01 


0.25 ± 0.03 


2.7 


0.37 


3.1 


r 

b 


p Oph 


147933 


194.5 ± 2.8 


58.9 ± 1.7 


16.7 ± 1.3 


23.5 ± 1.0 


62.9 ± 1.8 


17.9 


12.3 


5.8 


5.7 


91.6 ± 3.2 


98.3 ± 4.2 


0.30 ± 0.01 


0.27 ± 0.02 


1.5 


0.46 


2.4 


r 

b 


p Oph 


147955 


118.2 ±3.0 


62.8 ± 1.5 


10.6 ± 1.3 


22.6 ± 1.2 


35.0 ± 1.8 


12.2 


27. 1 


- 


- 


58.7 ± 1.7 


67.8 ± 3.4 


0.53 ± 0.02 


0.30 ± 0.04 


0.45 


0.49 


4.7 


b 


p Oph 


148041 


148.1 ±6.1 


50.6 ± 2.8 


25.9 ± 3.6 


48.0 ± 3.4 


94.4 ±3.1 


10.5 


20.0 


- 


5.2 


125.9 ± 6.3 


118.0 ±6.4 


34 + 02 


0.27 ± 0.04 


0.53 


0.40 


3.6 


r 

b 




148184 


81.7 ± 1.6 


57.7 ± 1.2 


14.2 ± 0.8 


19.0 ± 1.7 


45.5 ± 1.3 


23.6 


10.3 


2.0 


1.5 


72.9 ± 1.5 


85.1 ±3.2 


0.71 ± 0.02 


0.31 ± 0.02 


2.3 


0.72 


1.4 


£, 




148499 


105.4 ± 6.7 


34.1 ±4.0 


16.5 ±4.5 


24.7 ± 3.6 


25.8 ± 4.9 


6.4 


23.9 


- 


- 


104.0 ± 6.9 


110.9 ± 7.7 


0.32 ±0.04 


0.64 ± 0.21 


0.27 


0.36 


15.7 


<T 




148579 


104.6 ± 3.2 


55.0 ± 1.8 


13.8 ± 1.6 


18.4 ± 1.6 


39.9 ± 1.9 


11.1 


27.2 


- 


- 


54.3 ± 2.5 


65.0 ± 2.8 


0.53 ± 0.02 


0.35 ± 0.04 


0.41 


0.49 


5.8 


<T 


pOph 


148594 


78.8 ±3.0 


30.7 ± 1.9 


7.3 ± 1.1 


8.3 ± 1.1 


21.9 ± 2.1 


1.8 


7.8 




1.9 


39.6 ± 3.2 


54.6 ± 3.5 


0.39 ± 0.03 


0.33 ± 0.06 


0.23 


0.17 


10.2 






148605 


29.0 ± 1.2 


12.2 ± 1.6 


2.2 ±0.7 


< 1 


8.3 ± 1.1 


0.6 


1.2 






6.1 ± 1.6 


5.3 ± 1.8 


0.42 ± 0.06 


0.27 ±0.10 


0.50 


0.15 


7.4 




pOph 


148860 


108.9 ± 3.3 


28.7 ± 2.8 




15.3 ± 1.9 




2.0 






14.0 


62.3 ± 4.0 


78.4 ± 4.6 


0.26 ± 0.03 






0.14 


9.0 


ml 




149367 


108.2 ± 3.9 




6.6 ± 2.2 


9.7 ± 2.2 


33.4 ± 2.4 


4.2 


5.0 




3.1 


62.4 ± 4.9 


63.2 ±9.5 




0.20 ± 0.07 


0.84 


0.26 


4.8 






149387 


106.1 ±6.2 


32.0 ± 4.0 


16.4 ±3.1 


22.0 ± 3.0 




1 1.9 


1 1.2 




2.1 


90.7 ± 8.5 


95.9 ± 8.7 


0.30 ± 0.04 




1.1 


0.50 


9.8 


<r 




149438 


11.0 ±2.0 


15.5 ± 2.1 


<4 


<5.9 


< 5.9 










<4.2 




1.41 ±0.32 










•r 




149757 


51.0 ± 1.8 


41.0 ± 1.3 


7.7 ± 0.8 


18.9 ± 1.0 


45.7 ± 1.3 


18.7 


23.3 


6.2 




67.5 ± 2.1 


69.8 ± 1.9 


0.80 ± 0.04 


0.17 ±0.02 


0.8 


0.76 


1.8 


<r 




149883 


57.5 ± 3.0 


49.7 ± 3.0 


13.3 ±2.3 


33.6 ± 2.0 


33.0 ±2.1 


5.2 


10.0 




2.4 


65.5 ± 4.8 


81.5 ± 5.9 


0.86 ± 0.07 


0.40 ± 0.07 


0.52 


0.45 


5.3 


t 




149914 


116.0 ±2.3 


53.2 ± 1.3 


12.1 ±0.9 


26.9 ± 1.1 


57.0 ± 2.0 


18.1 


6.2 


10.2 


3.8 


109.9 ±2.1 


141.6 ±2.8 


0.46 ± 0.01 


0.21 ± 0.02 


2.9 


0.59 


2.5 


c 




150514 


114.2 ±3.6 


18.1 ± 2.5 


13.2 ±2.5 


20.5 ± 2.3 












45.2 ± 4.4 


55.2 ± 4.0 


0.16 ±0.02 










<r 




150814 


102.6 ± 4.5 


20.2 ± 3.0 


6.6 ± 2.0 


15.1 ±2.2 


27.2 ± 2.5 


2.9 


4.1 




1.8 


40.7 ± 3.6 


54.5 ± 6.4 


0.20 ± 0.03 


0.24 ± 0.08 


0.71 


0.20 


6.7 


o~ 




151012 


< 2.0 


<4.7 


<6.2 


<4.1 


< 11.0 










< 5.9 


















151346 


107.4 ±3.3 


55.8 ± 3.7 


12.1 ± 1.4 


26.7 ± 2.0 


69.8 ± 2.4 


26.9 


25.6 


8.7 




53.6 ±2.2 


58.9 ± 3.5 


0.52 ± 0.04 


0.17 ±0.02 


1.1 


0.69 


2.5 






151496 


148.1 ±4.8 


34.1 ± 2.7 


18.6 ±2.5 


32.3 ± 3.5 




3.3 


12.0 






57.8 ± 6.0 


65.8 ± 7.3 


0.23 ± 0.02 




0.28 


0.17 


9.9 


(T 




152516 


121.7 ±4.6 


26.6 ± 2.8 


12.4 ±2.1 


29.2 ± 2.7 


41.8 ± 2.6 










40.4 ±4.8 


53.6 ±7.0 


0.22 ± 0.02 


0.30 ± 0.05 








(T 




152655 


95.2 ± 4.0 


34.1 ± 2.3 


15.4 ±2.6 


22.6 ± 2.7 


44.4 ± 3.6 


3.3 


14.2 




1.7 


55.1 ± 5.5 


59.4 ± 5.7 


0.36 ± 0.03 


0.35 ± 0.07 


0.23 


0.24 


7.0 


( 




152657 


75.2 ±3.5 


29.0 ±3.1 


8.2 ± 1.4 


<4.9 


16.6 ± 2.1 


2.0 






2.9 


30.5 ± 3.7 




0.39 ± 0.04 


0.49 ±0.10 




0.19 


5.1 


( 





n 

Q 



C/3 

a 
a 



to 
— 



D 



—. i 



3 

c 

-3 

•a 



o 
-3 
c 

v. 

O 



5' 



■a 
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Table 5. Central radial velocities (in the heliocentric frame) for selected atomic and molecular lines. With an instrumental profile 
FWHM of about 6-7 km s _1 and pixels of 2 km s _1 , the accuracy of the measured radial velocity is 1-2 km s _1 . 



HD 


VCH+ 


VCH 


Radial velocities (km s" 

VCJV Veal v Kl 


- 1 ) 

VNal 




138503 


-6 


-6 


-6 


-6 


-6 




139094 


-8 


-9 


-11 -9 -9 


-8 


-10 




139160 


-7 


-7 


-8 


-8 


-9 




139486 






-13 


-13 


-13 




139518 










-15 




140543 


-10 




-13 


-13 


-14 




141180 










-7 


-17, +11 


141444 






-4 


-4 


-5 




141774 


-10 


-9 


-10 


-10 


-9 


-29, -18, +16 


142096 






-11 


-10 


-11 




142114 




-8 


-9 


-8 


-9 


-14 


142165 






-10 


-9 


-10 




142184 




-12 


-11 


-12 


-14 


-33 


142250 










-12 




142301 










-15 




142378 


-14 




-13 


-13 


-12 




142669 










-9 


-21 


142883 


-15 


-14 


-14 


-13 


-13 




142884 


-19 


-18 


-14 


-13 


-15 




142983 






-12 


-14 


-12 


-24 


142990 


-17 


-15 


-13 


-14 


-15 




143018 










-15 




143275 


-13 




-12 


-11 


-12 




143567 






-10 


-10 


-10 


-22 


143600 






-12 


-13 


-12 


-23 


143956 






-8 


-8 


-10 


-22 


144175 










-11 


-21 


144217 


-11 


-10 


-10 -10 


-10 


-10 


-23 


144218 


-11 




-10 -10 


-9 


-10 


-23 


144334 










-9 


-22 


144470 


-11 


-10 


-10 


-10 


-9 


-23 


144569 


-12 




-12 


-11 


-12 


-27 


144586 


-11 


-9 


-10 


-10 


-10 


-26 


144661 










-9 


-22 


144708 






-13 


-13 


-12 




144844 






-10 


-10 


-9 


-23 


144987 








-10 


-10 


-22 


145353 


-2 


-2 


-7 -4 


-5 


-5 


-21 


145482 








-1 1 


-1 


-20 


145483 








-11 




-25 


145502 


-11 


-9 


-10 -10 


-10 


-10 


-25 


145554 


-11 




-13 -10 


-10 


-10 


-25 


145556 


-3 


-4 


-4 


-4 


-4 


-25 
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Table 5. continued. 



Radial heliocentric velocities (km s 1 ) 



HD \'CH* VCH VCN v CaI Vki V^al 



145631 


-8 


-8 




-9 


-9 


-9 


-9 


-25 


145657 


-6 


-4 




-3 


-3 


-4 


-5 


-19 


145792 


-5 








-7 


-7 


-7 


-22 


145964 










-9 


-10 


-9 


-22 


146001 


-5 








-8? 


-9? 


-6 


-21 


146029 










-9 


-9 


-8 


-23 


146284 


-8 


-8 




-8 


-9 


-8 


-8 


-24 


146285 


-6 


-5 






-5 


-6 


-6 


-22 


146331 


-7 


-7 




-9 


-8 


-8 


-9 




146332 


-1 


-3 


-4 


-4 


-4 


-3 


-5 


+ 14 


146416 










-11 


-10 


-9 


-25 


146706 


-8 


-7 




-7 


-8 


-8 


-8 




147009 


-9 


-9 






-9 


-9 


-9 




147010 


-10 


-10 






-11 


-9 


-9 




147103 


-8 


-8 


-8 




-9 


-8 


-9 


-24 


147165 


-7 


-7 






-7 


-7 


-6 




147196 


-7 


-7 


-7 




-7 


-8 


-9 




147648 


-7 


-7 






-7 


-7 


-7 




147683 


-1 





-1 


-1 


-1 


-1 


+ 1 


-18 


147701 


-7 


-6 


-6 


-6 


-7 


-7 






147888 


-9 


-8 


-9 


-9 


-9 


-8 


-9 




147889 


-8 


-7 


-7 


-8 


-8 


-8 


-9 




147932 


-7 


-7 


-8 


-8 


-9 


-8 


-9 




147933 


-8 


-8 


-9 


-9 


-9 


-8 


-9 




147955 


-6 


-5 






-7 


-6 


-6 




148041 


-1 


-1 




-4 


-1 





-1 




148184 


-11 


-11 


-11 


-12 


-11 


-12 


-12 




148499 


-4 


-4 






-4 


-4 


-4 


+ 12 


148579 


-5 


-5 






-6 


-5 


-6 




148594 


-5 


-4 




-6 


-5 


-5 


-5 




148605 


-8 






-8 


-9 




-8 




148860 




-12 




-7 


-12 


-12 


-10 


-55 


149367 


-4 


-4 






-4 


-4 


-5 




149387 


-3 


-3 




-2 


-4 


-3 


-4 




149438 












-8? ? 


-8 




149757 


-15 


-15 


-15 




-15 


-15 


-15 


-27 


149883 


-5 


-4 




-6 


-4 


-4 


-5 




149914 


-12 


-11 


-12 


-10 


-11 


-11 


-10 




150514 


-8 








-7 


-7 


-9 


-27, +7 


150814 


-12 








-14 


-14 


-13 


+2 


151012 




-20? 






-19? 


-22? 






151346 


-7 


-7 


-7 




-7 


-7 


-8 




151496 


-9 


-8 






-9 


-9 


-10 




152516 


-7 








-6 


-6 


-7 




152655 


-10 


-9 




-8 


-10 


-10 


-9 




152657 














-4 





